Editora chefe 
Profº? Dr? Antonella Carvalho de Oliveira 
Editora executiva 
Natalia Oliveira 
Assistente editorial 
Flávia Roberta Barão 
Bibliotecária 
Janaina Ramos 
Projeto gráfico 
Bruno Oliveira 
Camila Alves de Cremo 
Daphynny Pamplona 2022 by Atena Editora 
Luiza Alves Batista Copyright O Atena Editora 
Natália Sandrini de Azevedo Copyright do texto O 2022 Os autores 
Imagens da capa Copyright da edição O 2022 Atena Editora 
iStock Direitos para esta edição cedidos à Atena 
Edição de arte Editora pelos autores. 
Luiza Alves Batista Open access publication by Atena Editora 


Todo o conteúdo deste livro está licenciado sob uma Licença de Atribuição 
Creative Commons. Atripuição-Não-Comercial-NãoDerivativos 4.0 


Internacional (CC BY-NC-ND 4.0). 


O conteúdo dos artigos e seus dados em sua forma, correção e confiabilidade são de responsabilidade 
exclusiva dos autores, inclusive não representam necessariamente a posição oficial da Atena Editora. 
Permitido o download da obra e o compartilhamento desde que sejam atribuídos créditos aos autores, 
mas sem a possibilidade de alterá-la de nenhuma forma ou utilizá-la para fins comerciais. 


Todos os manuscritos foram previamente submetidos à avaliação cega pelos pares, membros do 
Conselho Editorial desta Editora, tendo sido aprovados para a publicação com base em critérios de 
neutralidade e imparcialidade acadêmica. 


A Atena Editora é comprometida em garantir a integridade editorial em todas as etapas do processo 
de publicação, evitando plágio, dados ou resultados fraudulentos e impedindo que interesses 
financeiros comprometam os padrões éticos da publicação. Situações suspeitas de má conduta 
científica serão investigadas sob o mais alto padrão de rigor acadêmico e ético. 


Conselho Editorial 

Ciências Exatas e da Terra e Engenharias 

Prof. Dr. Adélio Alcino Sampaio Castro Machado - Universidade do Porto 

Profº Drº? Alana Maria Cerqueira de Oliveira — Instituto Federal do Acre 

Prof? Drº? Ana Grasielle Dionísio Corrêa - Universidade Presbiteriana Mackenzie 
Profº Drº? Ana Paula Florêncio Aires - Universidade de Trás-os-Montes e Alto Douro 
Prof. Dr. Carlos Eduardo Sanches de Andrade - Universidade Federal de Goiás 
Profº Dr? Carmen Lúcia Voigt - Universidade Norte do Paraná 


no 2022 


om [Atena 
MAR 


Prof. Dr. Cleiseano Emanuel da Silva Paniagua - Instituto Federal de Educação, Ciência e Tecnologia de Goiás 
Prof. Dr. Douglas Gonçalves da Silva - Universidade Estadual do Sudoeste da Bahia 

Prof. Dr. Eloi Rufato Junior - Universidade Tecnológica Federal do Paraná 

Prof? Dr? Érica de Melo Azevedo - Instituto Federal do Rio de Janeiro 

Prof. Dr. Fabrício Menezes Ramos - Instituto Federal do Pará 

Profº Dra. Jéssica Verger Nardeli - Universidade Estadual Paulista Júlio de Mesquita Filho 

Prof. Dr. Juliano Bitencourt Campos - Universidade do Extremo Sul Catarinense 

Prof. Dr. Juliano Carlo Rufino de Freitas - Universidade Federal de Campina Grande 

Profº Dr? Luciana do Nascimento Mendes - Instituto Federal de Educação, Ciência e Tecnologia do Rio Grande 
do Norte 

Prof. Dr. Marcelo Marques - Universidade Estadual de Maringá 

Prof. Dr. Marco Aurélio Kistemann Junior - Universidade Federal de Juiz de Fora 

Prof. Dr. Miguel Adriano Inácio - Instituto Nacional de Pesquisas Espaciais 

Profº Dr? Neiva Maria de Almeida - Universidade Federal da Paraíba 

Profº Drº? Natiéli Piovesan - Instituto Federal do Rio Grande do Norte 

Prof? Drº Priscila Tessmer Scaglioni - Universidade Federal de Pelotas 

Prof. Dr. Sidney Gonçalo de Lima - Universidade Federal do Piauí 

Prof. Dr. Takeshy Tachizawa - Faculdade de Campo Limpo Paulista 


- , [Ate 


na 


Ecdiirera 


Ano 2022 


Engenharia química: desenvolvimento de novos processos e produtos 2 


Diagramação: Camila Alves de Cremo 
Correção: Mariane Aparecida Freitas 
Indexação: Amanda Kelly da Costa Veiga 
Revisão: Os autores 
Organizador: Cleiseano Emanuel da Silva Paniagua 


Dados Internacionais de Catalogação na Publicação (CIP) 


EST Engenharia química: desenvolvimento de novos processos e 
produtos 2 / Organizador Cleiseano Emanuel da Silva 
Paniagua. - Ponta Grossa - PR: Atena, 2022. 


Formato: PDF 

Requisitos de sistema: Adobe Acrobat Reader 
Modo de acesso: World Wide Web 

Inclui bibliografia 

ISBN 978-65-258-0411-8 

DOI: https://doi.org/10.22533/at.ed.118221208 


1. Engenharia química. |. Paniagua, Cleiseano Emanuel 
da Silva (Organizador). Il. Título. 


CDD 660 
Elaborado por Bibliotecária Janaina Ramos - CRB-8/9166 


Atena Editora 

Ponta Grossa - Paraná - Brasil 
Telefone: +55 (42) 3323-5493 
www.atenaeditora.com.br 
contatodatenaeditora.com.br 


Atena 


Ebiicara 


áno 2022 


DECLARAÇÃO DOS AUTORES 


Os autores desta obra: 1. Atestam não possuir qualquer interesse comercial que constitua um conflito 
de interesses em relação ao artigo científico publicado; 2. Declaram que participaram ativamente da 
construção dos respectivos manuscritos, preferencialmente na: a) Concepção do estudo, e/ou 
aquisição de dados, e/ou análise e interpretação de dados; b) Elaboração do artigo ou revisão com 
vistas a tornar o material intelectualmente relevante; c) Aprovação final do manuscrito para 
submissão.; 3. Certificam que os artigos científicos publicados estão completamente isentos de dados 
e/ou resultados fraudulentos; 4. Confirmam a citação e a referência correta de todos os dados e de 
interpretações de dados de outras pesquisas; 5. Reconhecem terem informado todas as fontes de 
financiamento recebidas para a consecução da pesquisa; 6. Autorizam a edição da obra, que incluem 
os registros de ficha catalográfica, ISBN, DOI e demais indexadores, projeto visual e criação de capa, 
diagramação de miolo, assim como lançamento e divulgação da mesma conforme critérios da Atena 


Editora. 


Atena 


Ano 2027 


DECLARAÇÃO DA EDITORA 


A Atena Editora declara, para os devidos fins de direito, que: 1. A presente publicação constitui apenas 
transferência temporária dos direitos autorais, direito sobre a publicação, inclusive não constitui 
responsabilidade solidária na criação dos manuscritos publicados, nos termos previstos na Lei sobre 
direitos autorais (Lei 9610/98), no art. 184 do Código Penal e no art. 927 do Código Civil; 2. Autoriza 
e incentiva os autores a assinarem contratos com repositórios institucionais, com fins exclusivos de 
divulgação da obra, desde que com o devido reconhecimento de autoria e edição e sem qualquer 
finalidade comercial; 3. Todos os e-book são open access, desta forma não os comercializa em seu 
site, sites parceiros, plataformas de e-commerce, ou qualquer outro meio virtual ou físico, portanto, 
está isenta de repasses de direitos autorais aos autores; 4. Todos os membros do conselho editorial 
são doutores e vinculados a instituições de ensino superior públicas, conforme recomendação da 
CAPES para obtenção do Qualis livro; 5. Não cede, comercializa ou autoriza a utilização dos nomes e 
e-mails dos autores, bem como nenhum outro dado dos mesmos, para qualquer finalidade que não o 


escopo da divulgação desta obra. 


Atena 


Ano 2027 


APRESENTAÇÃO 


O e-book intitulado: “Engenharia química: Desenvolvimento de novos processos e 
produtos 2” é constituído por oito capítulos de livros que se distribuíram em três eixos- 
temáticos: í) produção e desenvolvimento de produtos de maior eficiência e de menor custo; 
ii) utilização de biomassa e controle de variável em diferentes processos industriais e; ill) 
geração de energia e estimativa de custo para recuperar espécies em recursos hídricos. 

O primeiro capítulo investigou a relação de um conjunto de filtros, utilizados em 
ambientes fechados, em relação à eficiência e o consumo energético resultando em um 
sistema que requeria um consumo de 8 W e uma eficiência de retenção de partículas 
abaixo do recomendado pela OMS. O capítulo 2 avaliou o desenvolvimento de produção 
de nanofibras via electrospinning a partir da análise de números adimensionais com 
múltiplas variáveis, os resultados sugerem estudos promissores que definiram a operação 
de produção de fibras pela via electrospinning. 

Os capítulos de 3 a 5 avaliaram: à) a importância do controle de temperatura no 
processo de esterilização de alimentos a partir de um modelamento matemático por meio 
de simulação computacional; ij) influência da temperatura na produção de ácido cítrico 
utilizando fungos da espécie Aspergillus Awamori e; iii) extração de xilose e glicose a partir 
do eucalipto (Eucalyptus sp.) a partir do pré-tratamento utilizando o ácido peracético e 
hidróxido de sódio. 

Por fim, os capítulos de 6 a 8 apresentaram trabalhos que avaliaram: i) utilização 
da energia térmica a partir de resíduos de Chillers, gerados em unidades hospitalares; 
ii) análise de custo para recuperar microalgas a partir de processos de microfiltração e; 
ili) remoção de arsênio (Ill) em matrizes aquosa empregando carbono pirolisado como 
adsorvente. 

Nesta perspectiva, a Atena Editora vem trabalhando de forma a estimular e incentivar 
cada vez mais pesquisadores do Brasil e de outros países a publicarem seus trabalhos 
com garantia de qualidade e excelência em forma de livros, capítulos de livros e artigos 


científicos 


Cleiseano Emanuel da Silva Paniagua 
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RESUMO: Uma das maiores preocupações 
sanitárias na atualidade é a contaminação por 
inalação de partículas ultrafinas. A solução para 
a mitigação deste risco encontra-se em parte 
nos sistemas de ventilação e condicionamento 
de ar (HVAC, de Heating, Ventilating and 
Air Conditioning). No presente trabalho, foi 
realizada uma análise dos parâmetros de projeto 
de um sistema HVAC composto por meios 
fitrantes de nanofibras de álcool polivinílico 
(PVA). Considerou-se uma sala contendo uma 
fonte de dispersão de gotículas de saliva com 
nanopartículas virais (SARS-COV-2). A sala 
também recebia partículas MP,, (Material 
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Particulado com diâmetro inferior a 2,5 um) do 
ambiente externo. Além dos filtros principais, 
foi considerado um filtro auxiliar com base no 
conceito de Minimum Efficiency Reporting Values 
(MERV). Foi utilizada como variável chave de 
análise a energia necessária para a ventilação 
da sala através da potência de ventiladores. 
Assim, foi possível atribuir uma estimativa de 
qual conjunto de filtros seria mais econômico em 
termos energéticos, visto que meios filtrantes de 
nanofibras não apenas possuem alta eficiência 
de coleta, como também alta queda de pressão. 
Verificou-se que o sistema mais econômico 
atingiu uma concentração de partículas abaixo 
do limite recomendado pela Organização Mundial 
de Saúde (OMS) após 167 minutos de operação. 
Este sistema requereria uma potência de 8 W 
considerando-se apenas a queda de pressão dos 
meios filtrantes limpos. Utilizando a metodologia 
empregada, foi possível obter um sistema que 
aliou relativo baixo consumo energético com 
a obtenção da concentração de aerossol em 
ambientes internos recomendada pela OMS. 
PALAVRAS-CHAVE: HVAC. Nanopartículas. 
Filtração de Ar. Electrospinning. Nanofibras 


DESIGN OF A VENTILATING AND 
AIR CONDITIONING SYSTEM USING 
NANOFIBROUS FILTER MEDIA 
ABSTRACT: One cf the biggest health concerns 
at present is the contamination with ultrafine 
particles through inhalation. The solution for 
the mitigation of this hazard is partially found in 
HVAC (Heating, Ventilating and Air Conditioning) 
systems. In this work, design parameters of an 
HVAC system composed of filter media made of 
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polyvinyl alcohol (PVA) nanofib rs were analyzed. It was considered a room with a source of 
saliva droplets containing viral nanoparticles (SARS-COV-2). The room also received outdoor 
PM, (Particulate Matter with diameter less than 2,5 um). Besides the main filters, an auxiliary 
filter was considered based on the concept of Minimum Efficiency Reporting Values (MERV). 
It was used as the key variable for the analysis the energy required for the ventilation of the 
room through the power consumption of fans. Therefore, it was possible to estimate which set 
of filters was the most cost-effective considering energy saving, since nanofibrous filter media 
present not only high collection efficienc , but also high pressure drop. It was verified that the 
most cost-effective system provided the particle concentration standard recommended by the 
World Health Organization (WHO) after 167 minutes of operation. This system would require 
8 W of power consumption considering only the pressure drop of the clean filter media. Using 
the presented methodology, it was possible to obtain a system that combined relatively low 
energy consumption with the obtention of an indoor particle concentration recommended by 
WHO. 

KEYWORDS: HVAC. Nanoparticles. Air Filtration. Electrospinning. Nanofibers 


11 INTRODUÇÃO 


Uma das aplicações imediatas dos meios filtrantes para filtração de ar é o emprego 
em sistema HVAC (Heating, Ventilation and Air Conditioning), que consiste em um conjunto 
de dutos para ventilação. Geralmente, tal sistema está relacionado ao conceito de Salas 
Limpas, no qual diferentes categorias de salas determinam a concentração máxima de 
particulado de determinado diâmetro no interior do ambiente, conforme ISO 146441-1 (TAN, 
2014). Neste tipo de sistema, a coleta de particulado pode ser realizada por meios filtrante 
ou precipitadores eletrostáticos, pois ambos possuem alta eficiência para partículas de 
pequeno tamanho (PARK et al., 2011; KIM et al., 2013). No entanto, filtros HEPA (High 
Efficiency Particulate Arrestance) são preferíveis, pois precipitadores podem produzir 
gás ozônio e consequentemente compostos orgânicos voláteis (VOCs — Volatile Organic 
Compounds), como terpenos e aldeídos. (WARING et al., 2008). 

O projeto de sistemas HVAC tem sido realizado para o controle da emissão de 
partículas virais em ambientes fechados, como no trabalho de Azimi e Stephens (2013), 
no qual se considerou um escritório fictício com dispersão de vírus Influenza. Mais 
recentemente no trabalho de Zhang et al. (2021), foi considerado um ônibus com dispersão 
de vírus SARS-COV-2. Este último vírus, bastante investigado após o início da pandemia 
de COVID (COrona Virus Disease) iniciada em 2019 (WHO, 2020), possui formato esférico 
e algum polimorfismo, com diâmetros variando entre 60 a 140 nm e apresentando spikes 
com 9 a 12 nm de tamanho em toda a sua superfície externa (ZHU et al., 2019; SCHELLER 
et al., 2020). 

No geral, partículas patogênicas são consideradas como particulados poluentes de 
ar, visto que se suspendem em ar facilmente: vírus possuem tamanhos entre 20 e 400 
nm, bactérias, entre 200 nm e 2,0 um, enquanto esporos de fungos possuem entre 2,0 
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e 8,0 um (JOHNSON et al. 2009; SUNG et al., 2016; ZHOU et al., 2018). Uma forma 
comum de dispersão destes agentes é através da emissão de gotículas a partir de um 
indivíduo contaminado, ao tossir, espirrar ou falar, emitindo gotículas entre 12 e 21 um 
contendo enzimas, eletrólitos e material genético destes agentes (STADNYTSKYI et al., 
2020). No entanto, a dinâmica de aerossóis é complexa, visto que há possibilidade de 
partículas virais dispersas coalescerem até tamanhos maiores que 2,5 um (SETTI et al., 
2020). Adicionalmente, é possível que gotículas de saliva como as que contém SARS- 
COV-2 reduzam seu tamanho após dispersas em ar devido à evaporação da água, o que é 
dependente da temperatura e umidade ambiente (ROHIT et al., 2020). 

Portanto, o presente estudo consistiu no projeto de sistema HVAC para a retenção 
de MP,, em uma sala contendo uma fonte de dispersão de gotículas de saliva com 
partículas virais (SARS-COV-2) e considerando-se um filtro auxiliar com base no conceito 
de Minimum Efficiency Reporting Values (MERV) (ASHRAE, 2007; STEPHENS; SIEGEL; 
2012). Foi utilizada como variável final de análise a energia necessária para a ventilação 
do sistema através da potência de ventiladores. Assim, foi possível estimar qual conjunto 
de filtros seria mais econômico em termos energéticos, uma vez que a presença de um 
meio filtrante de nanofibras não apenas aumenta a eficiência de coleta, como também a 
queda de pressão da operação. Foi fixada a concentração de saída de MP, , da sala como 
a concentração limite média diária estipulada pela Organização Mundial de Saúde (OMS), 
igual a 15 ug/mº para ambientes internos (indoor) (2021). Este valor limite foi recentemente 
atualizado, o que implica na exigência de sistemas de purificação de ar mais eficientes. Isto 
ressalta a importância do presente trabalho no que tange a demonstrar a aplicabilidade de 
sistemas de filtração de ar de alta eficiênci 


21 METODOLOGIA 


A Figura 1 apresenta uma ilustração do balanço de massa para um sistema HVAC 
(Box model), contendo um reciclo e matematicamente expresso por (TAN, 2014): 


Px 


= [Oo x Co xl ER+ Gem Cx (1 ER m]- O, x E + 0x0 x Ef) (1) 


E fr 


em que Vé o volume interno útil da sala de dimensões X, Ye Z. O parâmetro Q, é a 
vazão volumétrica de ar, cujo fluxo é representado pela seta com preenchimento em preto 
da Figura 1 (setas com preenchimentos com diferentes gradações de cinza são ilustrativas 
da redução da concentração de aerossol). Esta corrente de ar contém particulado com 
concentração C,, que atravessa um filtro externo (make up) com eficiência de filtraçã 
Ef,. A variável Q se refere à parte da vazão de saída Q, reciclada e que atravessa um 
filtro de reciclo de eficiência Ef, antes de retornar ao duto de ventilação (não explicitado 
na Figura 1 para simplificar o esquema). A variável m corresponde à taxa de geração de 


particulado gerada no interior da sala, aqui a taxa de geração de gotículas contendo carga 
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viral. O parâmetro C, se refere à concentração de particulado interna da sala (e de saída 
considerando o ar perfeitamente misturado), resultado da concentração de partículas antes 
do início da operação de ventilação e das vazões de particulado de entrada e de reciclo, 
sendo, portanto, dependente do tempo de operação t. Por fim, Ef, se refere à eficiência do 
meio filtrante no interior da sala, representado sob a forma de linhas tracejadas na Figura 


1, do mesmo modo como as correntes de entrada e saída de ar desse filtro 


“= da recirculação 


Figura 1. Ilustração de um sistema HVAC em uma sala 


Fonte: Acervo pessoal 


A solução do equacionamento resulta em: 
E = Cal x erpl=titi+ O. (2) 


, 
[e Lg 


Co x (L= ELA mido = 0 x l1=ERI+O= EL) [3] 


T= Fido = O x(l=EL)+0 = Ef) (4) 


sendo C, a concentração inicial de particulado no interior da sala, C, a concentração 
de partículas no interior da sala quando o sistema atinge o estado estacionário e rt, uma 
“constante de tempo”, ou seja, um valor característico do sistema com unidade de tempo e 
que permite caracterizar a concentração transiente no espaço (TAN, 2014). 

Para a execução do projeto, considerou-se uma sala contendo apenas um indivíduo, 
contaminado com SARS-COV-2 emitindo gotículas de saliva a uma taxa de 50 gotas por 
segundo, com mediana de -0,8 um de diâmetro e desvio padrão geométrico dependendo 
da amplitude da vocalização (ASADI et al., 2019). A carga viral de SARS-COV-2 na saliva 
de uma pessoa infectada na primeira semana após o aparecimento dos primeiros sintomas 
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varia entre 10º a 10º partículas/mL (ZHU et al., 2020). Considerando a concentração máxima 
de vírus e a densidade da saliva igual a 1000 kg/mº, foi calculada a taxa de emissão de 
vírus no interior das gotículas para cada tamanho de gota. Sendo a faixa de tamanho de 
partícula de SARS-COV-2 de 60 a 140 nm, foi estipulada uma distribuição com mediana 
igual a 100 nm e desvio padrão geométrico de 1,1. Tomando-se a densidade do vírus como 
igual a 1190 kg/m? (HUANG et al., 2006), foi estimada a taxa mássica de geração de vírus 
no interior da sala. Assim, m foi obtida a partir da soma das taxas de geração de gotas 
e de vírus a partir de um único indivíduo ocupante da sala, contaminado com o vírus e 
vocalizando a uma taxa máxima de emissão de gotículas, as quais continham carga viral 
máxima. 

O valor de Q necessário para ventilar a sala foi calculado com base na especificaçã 
da ASHRAE 62.1 (2010), igual a 8,5 m?/h por pessoa + 1,1 m*/h por área de piso, que foi 
estipulada como igual a 25 m?. Assim, Xe Y eram iguais a 5 me Z iguala 3 m. A área 
de filtração para cada filtro foi estipulada como igual 0,6x0,6 (m?), conforme realizado por 
Azimi e Stephens (2013). 

A Tabela 1 apresenta valores de classificação de filtros de acordo com o conceito 
de Minimum Efficiency Reporting Values (MERV) (ASHRAE, 2007; STEPHENS; SIEGEL; 
2012). 


MERV 0,34 um(%) 1-3 um (%) 3-10 um (%) APmédia: (Pa) 

4 1 9 15 73 

7 17 46 50 11 

1º 30 65 85 141 

13º 70 90 90 144 

14º 80 90 90 188 

15º 90 90 90 159 

16: 95 95 95 157 
HEPAº 99,9 99,9 99,9 374 


Tabela 1. Classificação de elementos filtrant 


Fonte: “ASHRAE, 2007; "STEPHENS; SIEGEL; 2012; cAZIMI, STEPHENS, 2013 


Desconsiderando-se a presença de nanopartículas virais advindas do exterior, 
estipulou-se que o filtro make up seria constituído de um filtro genérico que atendesse 
às classificações mostradas na Tabela 1. A Tabela 1 também apresenta valores de queda 
de pressão média obtidos por Azimi e Stephens (2013) para meios filtrantes comerciais 
enquadrados nestas classificações. Este parâmetro foi utilizado para o cálculo da energia 
requerida para ventilação. 

Salienta-se que a denominação “filtro HEPA” a ser utilizada nesta análise corresponde 


à denominação adotada no trabalho de Azimi e Stephens (2013). Na realidade, de acordo 
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com a EN 1822, meios filtrantes com esta eficiência são classificados como E12, enquanto 
a ISO 29463 os classifica como ISO 30E (EMW FI TERTECHNIK, 2022). 

Os filtros interno e de reciclo, em contato com as gotículas e partículas virais, seriam 
compostos dos meios filtrantes de nanofibras. Foram utilizados resultados experimentais de 
caracterização morfológica e de performance de meios filtrantes de nanofibras produzidos 
por electrospinning, obtidos em trabalho prévio (OLIVEIRA et al., 2021). No estudo, foram 
investigadas diferentes composições de soluções poliméricas de álcool polivinílico (PVA) 
alterando-se a concentração do surfactante Triton X-100, utilizado para se obter fibra 
nanométricas para filtração de nanopartículas. 

É preciso ressaltar que, para efeitos de simplificação, o valor de queda de pressão 
tomado para os meios filtrantes de PVA foi o de queda de pressão inicial (de filtro limpo), 
visto não haver dados da performance desses meios filtrantes após longos períodos de 
filtração. O cálculo da energia requerida foi bastante simplificado, considerando a eficiênci 
dos ventiladores e dos motores iguais a 1 (um) apenas para efeitos comparativos. A potência 
é expressa por (AZIMI; STEPHENS; 2013): 


Pot = (0, 4x AP + O, E AP aaa (5) 


Nr] , 


Para a análise do gasto energético (queda de pressão) correspondente aos meios 
fitrantes de nanofibras, como função do tempo de fiação e da concentração de surfactante 
utilizado na solução polimérica, foi utilizada a definição de g amatura (W): 


LATER [Er 


sendo a e L respectivamente a densidade de empacotamento (contraparte da 


porosidade) e a espessura do meio filtrante de nanofibras; , é a densidade do polímero. 


Considerando-se a área de filtração utilizada nos féses experimentais, que possuía 
seção circular de diâmetro interno igual a 1 polegada (OLIVEIRA etal., 2021), e considerando- 
se homogênea a deposição espacial das camadas de fibras por electrospinning, obteve-se 
a massa total de fibras depositada para aquela área útil de meio filtrante. Considerando-se 
uma relação linear entre a massa depositada e o tempo de fiação, obteve-se a taxa de 
deposição de nanofibras para cada composição de solução considerando-se o tempo de 
fiação experimental, igual a 15 min (OLIVEIRA et al., 2021). A Tabela 2 exibe os valores 


calculados. 
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suriactante Cemlm) fbra(mm) 90) (tamin)  Massa(60) — goposição (kg) 
0 234 0,047 11,6 3,59E-07 4,0E-10 
0,25 229 0,047 12,8 3,96E-07 4,4E-10 
0,50 186 0,043 11,7 3,31E-07 3,7E-10 
0,75 191 0,052 10,0 3,43E-07 3,8E-10 
1,00 180 0,054 9,38 3,34E-07 3,7E-10 


Tabela 2. Caracterização dos elementos filtrantes de Oliveira et al. (2021) e taxas calculadas de 
deposição de fibra 


A partir de tais taxas (Tabela 2), foi possível obter a espessura de meio filtrante como 
função do tempo de fiação para cada tipo de elemento filtrante produzido, bem como obter 
a queda de pressão a partir da equação de Davies (1973), considerando que a porosidade/ 
densidade de empacotamento independe do tempo de fiação. Esta e uação é dada por: 


dP=to E px loeló (+ ha!) (7] 


na qual D é o diâmetro da fibra e ve Ho são respectivamente a velocidade superficia 
e a viscosidade do ar. 

Usando-se a velocidade de filtração utilizada nos experimentos, igual a 5,0 cm/s 
(OLIVEIRA et al., 2021), Q, foi obtida como igual a 0,05x0,6x0,6=0,018 m*/s. Tal velocidade 
foi utilizada uma vez que os dados de performance de filtração daquele trabalho (-95% de 
eficiência global para filtração de partículas nanométricas) foram obtidos nesta condição. 
Estipulou-se que Q=Q, /3, conforme realizado por Tan (2014). As concentrações de MP, , 
do ambiente externo (outdoor) (C) e interno anteriores ao início da ventilação (C,) foram 
tomadas como iguais a 155 e 118 ug/mº, relativas a um apartamento em perímetro urbano 
(DENG et al., 2017). Foi então estipulada uma distribuição de MP, , outdoor com mediana 
de 1,4 um e desvio de 1,1. Como anteriormente mencionado, como a concentração limite 
média diária estipulada pela OMS é igual a 15 ug/mº* para ambientes indoor (2021), 
estipulou-se que este deveria ser o valor limite de C,., abaixo do qual o sistema HVAC seria 
eficaz para a ventilação da sala em termos de concentração de p rticulado. 


31 RESULTADOS E DISCUSSÃO 


A Figura 2 apresenta os dados de função de densidade de probabilidade (ver TAN, 
2014) do particulado calculada para cada tipo de meio filtrante candidato a filtro de make 
up, considerando-se um particulado advindo do ambiente externo MP,, (mediana: 1,4 
um; desvio padrão geométrico: 1,1). À esquerda (a), para o aerossol na entrada do meio 
filtrante e à direita (b), na saída, para cada tipo de filtro MERV e para filtro HEPA. Uma 
análise rápida permite estabelecer relações com resultados posteriores. 
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Figura 2. Função densidade de probabilidade do particulado antes (a) e após (b) a filtração com 
diferentes classificações para o filtr make up 


Verifica-se não apenas a redução da densidade de probabilidade do aerossol com 
o aumento da classificação MERV devido ao aumento da eficiência (Tabela 1), como 
também a drástica redução obtida entre os meios filtrantes classificados por Stephens e 
Siegel (2012) e os da norma ASHRAE 52.2. A propósito, constata-se uma modificaçã 
na mediana da distribuição das curvas a depender da classificação do meio filtrante. Isto 
ocorreu porque, para os tipos MERV 15 e 16 e para o filtro HEPA, a eficiência de coleta 
é igual para as subdivisões da distribuição de tamanho de partícula (Tabela 1), o que não 
ocorre com os demais. 

A partir dos cálculos e hipóteses assumidas, obteve-se os valores de C, (Figura 3 
a)er (Figura 3 b) para os sistemas contendo diferentes filtros make up e filtros interno e 
de reciclo de eficiências Ef= Ef, conforme esquema da Figura 1. A concentração de saída 
estipulada (denominada limite) está destacada em vermelho para o gráfico de C,.. 
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Figura 3. Valores de C, (a) e 1 (b) para diferentes configurações de filt make up e filtros de reciclo e 
interno do sistema HVAC proposto 
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Da Figura 3 (a), verifica-se que, para as condições operacionais e parâmetros de 
entrada utilizados neste estudo, não há combinação viável entre meios filtrantes internos 
e de reciclo de quaisquer eficiências com filtros make up com classificação MERV 4, 7, 
1 e 13. Para a faixa de eficiência observada (50 a 99,99%), a concentração de MP, 
(incluindo gotículas e partículas virais) ao atingir-se o estado estacionário após o início da 
operação de ventilação é superior ao limite de 15 ug/m* recomendado pela OMS (2021). 
Para o filtro MERV 14, sua associação com os meios filtrantes internos e de reciclo passa 
a ser viável para eficiências de filtração destes acima de 70%. Contudo, é necessário frisar 
que a eficiência de filtração para nanopartículas é fundamental no caso de retenção de 
partículas virais em ambientes fechados. Um sistema de ventilação e condicionamento de 
ar inadequado pode simplesmente ajudar na dispersão do vírus em vez de retê-lo, como 
reportado em um restaurante em Guangzhou (China) contendo um cliente contaminado 
com SARS-COV-2 (LU et al., 2020). Assim, além de atender as normas de concentração 
da Organização Mundial de Saúde, é necessário um sistema eficiente na mitigação da 
dispersão destas partículas nanométricas. Desta forma, foram considerados neste estudo 
para compor os filtros interno e de reciclo os meios filtrantes de nanofibras de PVA 
desenvolvidos por Oliveira et al. (2021), uma vez que forneceram eficiências de coleta de 
nanopartículas (5,94-224,7 nm) de -95%. 

Verifica-se na Figura 3 (a) que o filtro HEPA utilizado como filtro make up faz ser 
dispensável o uso de meios filtrantes para atingir-se a concentração de saída estipulada. 
Nota-se que os resultados obtidos para C, estão bastante abaixo do limite delimitado e até 
mesmo próximos a zero, independentemente da eficiência dos filtros interno e de reciclo. 
Os meios filtrantes tipo MERV 15 e 16 possuem quase o mesmo comportamento, exceto 
que apresentam pequena variação de C,. com o aumento de Ef=Ef. Pondera-se que 
elementos filtrantes tão eficientes não sejam necessários neste sistema à primeira vista e 
até mesmo que encareçam o processo em termos de consumo energético. Isto pode ser 
verificado na Tabela 1 ao comparar-se a queda de pressão média do filtro tipo HEPA em 
relação aos demais, seguidamente pelos filtros tipo MERV 15 e 16. Ainda na Figura 3 (b), 
nota-se a redução da variável 7 com o aumento da eficiência dos meios filtrantes dos filtro 
interno e de reciclo e que não envolve o tipo de filtro de make up. 

A Figura 4 exibe a relação da concentração no ambiente interno com tempo de 
operação para as diferentes configurações possíveis de filtro make up e filtros interno 
e de reciclo, com diferentes eficiências. Como complemento dos resultados da Figura 
3, os gráficos da Figura 4 permitem visualizar quais configurações de filtr make up e 
filtros interno e de reciclo, a partir de determinado tempo de operação do sistema HVAC, 
atingem a concentração limite de MP,, recomendada pela OMS (2021) (em vermelho e 
denominada limite). Novamente, as configurações utilizando filtros de classificação MERV 
4,7,11 e 13 não atingem a concentração de saída estipulada mesmo após 10.000 s da 
partida da operação de ventilação (1,2-1,9 vezes a “constante de tempo” a depender de 
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Ef= Ef, conforme Figura 3). Configurações utilizando filtro tipo MERV 14 como filtro make 
up atingirão tal valor após 10.000 s para eficiências acima de 90%; demais configuraçõe 
atingirão tal limite de concentração a tempos cada vez menores, a depender da eficiênci 
dos meios filtrantes utilizados. Salienta-se que a faixa de tamanho de partícula analisada 
neste estudo é mais extensa que aquela avaliada no trabalho de Oliveira et al. (2021) 
no qual foram baseados os dados de performance dos meios filtrantes de PVA. Naquele 
estudo, foram atingidas eficiências de -95% para partículas que inclusive constituem a 
faixa de mínimo de eficiência de filtração (HINDS, 1998). Portanto, é esperado que os 
meios filtrantes de PVA utilizados neste projeto coletem com alta eficiência as partículas da 
faixa avaliada, ou seja, as partículas MP, , advindas do exterior, além das nanopartículas 
geradas no ambiente interno. 
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Figura 4. Valores de C, em função do tempo de operação para diferentes configurações d sistema HVAC 
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A decisão sobre a escolha do meio filtrante adequado deve ser dada, além da 
eficiência, pelo consumo energético requerido do processo, além de custos relacionados à 
aquisição e troca dos meios filtrantes. Estes últimos custos envolvem aspectos comerciais 
que estão fora do escopo deste trabalho. Assim, a Figura 5 apresenta os resultados de 
potência requerida para a ventilação do sistema HVAC ilustrado na Figura 1 com diferentes 
configurações de filtros make up (apenas para MERV 14, 15 e 16), além dos filtros interno 
e de reciclo. 

Nota-se na Figura 5 que, como esperado, o aumento do tempo de deposição dos 
meios filtrantes utilizados — e consequentemente sua eficiência — implica no aumento do 
consumo de energia requerida para a operação do sistema HVAC, o que está relacionado 
com a permeabilidade dos meios filtrantes e consequentemente à resistência que estes 
meios impõem à passagem da corrente gasosa (DAVIES, 1973). Considerando, por 
exemplo, um sistema HVAC com um meio filtrante MERV 15 como filtro make up e o restante 
dos filtros compostos por meios filtrantes fiados a partir de soluções de 0,25; 0,50; 0,75 e 
1,0% de PVA, que obtiveram experimentalmente -95% de eficiência para nanopartículas 
da faixa de tamanho que inclui a faixa de vírus tipo SARS-COV-2, a potência requerida para 
ventilação seria respectivamente de 9; 11; 11 e 12 W para tempos de fiação de 30 min, 
ao passo que a substituição por um meio filtrante tipo MERV 14 implicaria num consumo 
de 8; 9; 10 e 10,8 W mantidas as demais condições. Por outro lado, a maior quantidade 
de particulado a adentrar a sala partindo de um filtro make up com menor capacidade de 
retenção de MP, , implicaria em uma saturação mais rápida dos meios filtrantes dos filtro 
interno e de reciclo, que poderia exigir trocas mais frequentes destes. 

Além do consumo energético, custos associados à compra e troca dos elementos 
fitrantes devem ser considerados na escolha das configurações mais econômicas. Azimi 
e Stephans (2013), ao projetar o sistema HVAC para a retenção de partículas virais de 
Influenza partindo-se dos meios filtrantes da Tabela 1, realizaram tais considerações. No 
entanto, não é possível mensurar facilmente o custo de compra de elementos filtrante 
desenvolvidos em laboratório, sendo possível, em um primeiro momento, compará-los 
quanto a aspectos de performance. 
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Figura 5. Potência total requerida para ventilação dos sistemas HVAC possíveis (cada gráfico 
corresponde ao meio filtrante de P A produzido com distintos tempos de fiação para filtros interno e d 
reciclo — vertical — em associação com filtros ME V para filtros de make up 


Verifica-se, a priori que, para a retenção de -95% do particulado gerado no interior 
da sala e advindo do ambiente externo sob as condições avaliadas do processo, é mais 
economicamente viável o uso de um meio filtrante tipo MERV 14 no filtro make up e filtro 
interno e de reciclo contendo elementos filtrantes de PVA fabricados a partir de solução 
com 0,25% em massa de PVA. Nestas condições, o tempo necessário para obter um 
ambiente com concentração indoor de MP,, limite de acordo com a OMS seria de 167 
minutos após a partida da operação. Para tempos menores, ou seja, caso se deseje ocupar 
mais rapidamente a sala com a concentração de aerossol recomendada pela OMS, uma 
nova combinação teria que ser realizada com um filtro de make up de maior eficiência 
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mantidas as demais condições. 

Há ainda que ser considerada a retenção das partículas virais nos elementos 
filtrantes dos filtros. Apesar de não ser o foco inicial no desenvolvimento do material, a 
aplicação em sistemas de ventilação pode sugerir a futura aplicação de agentes biocidas 
na produção das nanofibras de PVA, assim como outros trabalhos da literatura têm se 
especializado na aplicação de nanopartículas nas nanofibras, como as de prata em meios 
filtrantes destinados à filtração de ar contaminado com agentes biológicos (ROSA et al., 
2017; BORTOLASSI et a!l., 2019; TREMILIOSI et a!., 2020). De toda forma, o cuidado com a 
limpeza e a validade de meios filtrantes, seja em sistemas de ventilação e condicionamento 
de ar de edifícios ou de automóveis, deve ser sempre tomado de forma a evitar o crescimento 
de agentes biológicos retidos nas superfícies das fibras de meios filtrantes, que ocorre 
dependendo das condições ambientais do processo (KELKAR et al., 2005; LI et al., 2013; 
GOLOFIT-SZYMCZAK et al., 2019). 


41 CONCLUSÕES 


Foi possível elaborar o projeto de um sistema HVAC para o controle da concentração 
indoor de particulado MP, em uma sala na qual nanopartículas virais de SARS-COV-2 
eram geradas. Para isto, foi projetado um conjunto de meios filtrantes comerciais com os 
meios filtrantes de nanofibras de PVA produzidos experimentalmente em trabalho prévio. 
Foram considerados aspectos energéticos (queda de pressão) envolvidos na operação 
com meios filtrantes comerciais de distintas eficiências e também os relativos aos meios 
filtrantes de PVA produzidos com diferentes concentrações de surfactante Triton X-100. Foi 
possível obter um sistema que aliasse relativo baixo consumo energético com a obtenção 
da concentração de aerossol indoor recomendada pela OMS. Por exemplo, um sistema 
composto de um filtro MERV 14 (filtro de make up) e meios filtrantes de PVA com 0,25% 
m/m de Triton X-100 (filtros interno e de reciclo) seria adequado para manter a sala sob 
a concentração recomendada após 167 min do início da operação. Espera-se que esse 
sistema retenha o particulado nanométrico gerado na sala graças à eficiência de coleta das 
nanofibras que foi comprovada experimentalmente. Este sistema requereria uma potência 


de 8 W considerando-se apenas a queda de pressão dos meios filt antes limpos. 
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RESUMO: Números  adimensionais são 
facilitadores para a descrição matemática de 
relações entre diferentes forças envolvidas em 
diversos fenômenos, bem como para a realização 
de experimentos e análises com múltiplas 
variáveis. No presente trabalho, são apresentados 
adimensionais desenvolvidos a partir do teorema 
de Vaschy-Buckingham, considerando o processo 
de produção de nanofibras via electrospinning. 
No desenvolvimento, as propriedades físico- 
químicas do fluido são colocadas em destaque. 
Os adimensionais obtidos foram utilizados 
com dados experimentais da literatura e suas 
possíveis aplicações foram exploradas em um 
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DE ELECTROSPINNING 


estudo teórico sobre filtração de nanopartículas. 
Foi possível utilizar o adimensional e os dados 
experimentais da literatura para predizer a 
performance de meios filrantes de álcool 
polivinílico (PVA) produzidos com diferentes 
tempos de fiação para diferentes distribuições 
de tamanho de partícula. Com a utilização dos 
adimensionais, futuros trabalhos poderão definir 
regimes de operação e nos quais poderá ser 
possível obter fibras a partir de electrospinning 
de acordo com as propriedades da solução, 
bem como definir regimes de ocorrência de 
nódulos (beads) nas fibras de acordo com estas 
propriedades. 

PALAVRAS-CHAVE: Números adimensionais. 
Filtração de Ar. Electrospinning. Nanofibras. 
Nanopartículas. 


OBTENTION OF DIMENSIONLESS 
NUMBERS FOR THE STUDY OF 
ELECTROSPINNING 


ABSTRACT: Dimensionless numbers are useful 
to describe relationships between different 
forces involved in several phenomena, besides 
their use for the execution of experiments 
and analyses with several variables. The 
dimensionless numbers presented in this work 
were developed using the Vaschy-Buckingham 
theorem considering the process of nanofibers 
production via electrospinning. In the theoretical 
development, the physical-chemical properties 
of the polymer fluid were highlighted. The 
dimensionless numbers obtained were used with 
experimental data from the scientific literature 
and their possible applications were explored in 
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a theoretical study on nanoparticle filtration. It was possible to use the dimensionless number 
and the experimental data to predict the performance of polyvinyl alcohol (PVA) filter media 
produced with different spinning times for different particle size distributions. By using the 
dimensionless numbers, further works may determine operation regimes in which the obtention 
of fibers through electrospinning are viable in accordance with the solution properties. They 
also may determine operation regimes in which beads appear in the fibers in accordance with 
these properties. 

KEYWORDS: Dimensionless Numbers. Air Filtration. Electrospinning. Nanofibers. 
Nanoparticles. 


11 INTRODUÇÃO 


A filtração de gases é um dos processos mais exitosos de separação gás-sólido 
para a mitigação da emissão de particulados à atmosfera, sendo possível obter eficiência 
de coleta próximas a 100% em uma ampla faixa granulométrica (SUTHERLAND; CHASE, 
2008). Sua operação baseia-se na atuação conjunta de uma série de mecanismos físicos 
que culminam na aderência da partícula dispersa no fluido sobre a fibra do meio filtrante 
Contudo, a coleta de nanopartículas é dificultada devido à sua massa desprezível (que 
resulta em baixas eficiências dos mecanismos gravitacionais e inerciais) e pequeno volume 
(baixa eficiência de interceptação direta) (HINDS, 1998) 

A confecção de meios filtrantes contendo nanofibras através do processo de 
electrospinning tem sido uma aposta tecnológica para contornar este problema. Este 
método consiste na passagem de fluido polimérico por um capilar conectado a uma fonte 
de alta tensão, que gera uma diferença de potencial entre o capilar e um coletor localizado 
a uma distância específica. As forças elétricas advindas da diferença de potencial, a partir 
de uma voltagem crítica, superam a tensão superficial do fluido. Isto forma um jato que se 
prolonga de uma formação cônica conhecida como cone de Taylor e que se estende até 
atingir o coletor, produzindo um emaranhado de fibras com a volatilização do solvente que 
solubilizava o polímero, no caso de o fluido polimérico ser uma solução (RENEKER; CHUN, 
1996). 

Como forma de abordagem do complexo processo de electrospinning de soluções 
poliméricas, o emprego de números adimensionais com as propriedades destes fluido 
(condutividade elétrica, tensão superficial e viscosidade) e variáveis de processo (vazão 
de alimentação, voltagem aplicada, distância entre agulha e coletor, rotação do coletor, 
campo elétrico, condições ambientais etc.) poderia reduzir a quantidade de variáveis a 
serem avaliadas. Na literatura científica, números adimensionais que foram obtidos a partir 
do teorema dos grupos 11 de Vaschy-Buckingham têm alcançado sucesso na minimização 
do número de variáveis avaliadas, simplificando os esforços na aquisição de resultados 
teóricos e experimentais para processos sob a influência de múltiplos efeitos (WELTY, 2008). 


Alguns trabalhos envolvendo electrospinning, de fato, produziram números adimensionais 
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para descrever especificamente a formação e instabilidade dos jatos de polímero e as 
características das fibras produzidas (FENG, 2001; McKINLEY, 2005; GADKARI, 2014; 
YARIN etal., 2014). Contudo, não foi encontrado um trabalho que sumarizasse as principais 
variáveis de processo e as propriedades do fluido, desde a formação do cone de Taylor até 
a evaporação do solvente durante a deposição das fibras no coletor, especialmente no caso 
de coletores rotativos. Salienta-se que este tipo de aparato promove uma fluidodinâmic 
distinta relacionada à deposição das fibras, ao seu alinhamento na trama de acordo com 
a velocidade de rotação e à maior taxa de evaporação de solvente quando se compara 
com processos usando-se coletores fixos (RAMAKRISHNA et al., 2005; MEDEIROS et al., 
2008). 

Neste trabalho, será apresentado um desenvolvimento de adimensionais a partir 
da teoria dos grupos Tr de Vaschy-Buckingham, que poderá ser demonstrado como útil 
em estudos envolvendo o processo de electrospinning, sendo este um processo que 
sofre influência de inúmeras variáveis de alimentação e operacionais. Os adimensionais 
desenvolvidos neste trabalho serão aplicados utilizando-se dados experimentais obtidos 
na literatura para a predição da eficiência de meios filtrantes de nanofibras para a filtraçã 
de nanopartículas. Com os resultados deste trabalho, espera-se que estudos posteriores 
possam averiguar demais aplicações dos adimensionais obtidos para que, em caso de 


êxito, facilitem futuros estudos sobre este processo. 


21 DESENVOLVIMENTO TEÓRICO 


Para a obtenção dos adimensionais, foi realizado um levantamento das principais 
variáveis que influenciam o processo de electrospinning quanto às condições de alimentação 
(propriedades da solução polimérica) e condições operacionais (vazão de alimentação de 
solução, voltagem, diâmetro da agulha/capilar etc) (RAMAKRISHNA et al., 2005; YARIN 
et al.; 2014). Assim, a Tabela 1 apresenta as variáveis consideradas para a aplicação do 
teorema Tr de Vaschy-Buckingham (WELTY, 2008) baseando-se nos resultados recorrentes 


da literatura. 


; az Dimensô istema Internacional 
Símbolo Variável ensões do Sistema Internacional de 


Unidades 
q Diâmetro da fibra IL] 
o Densidade do fluido IM L3] 
p” Viscosidade do fluido IM LT] 
Y Condutividade do fluido [A2 Tº Mº 13] 
fo Tensão superficial do fluido IM T2] 
e Voltagem IM LZ2 A? T3] 
E: Distância entre agulha e coletor IL] 
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n Diâmetro da agulha [L] 

q Vazão do fluido [1º 7] 
4 Rotação do coletor [LT] 
w Pressão de vapor do fluido IML“ T2 


Tabela 1. Variáveis consideradas para a obtenção dos grupos tr 


Da Tabela 1, verifica-se a existência de 11 variáveis principais (sendo a variável w' 
designada como dependente das demais) e 4 dimensões (comprimento L, massa M, tempo 
Te intensidade de corrente elétrica A). Portanto 7 grupos Tr possíveis. Aqui, não estão sendo 
consideradas propriedades adimensionais, como as permissividades elétricas relativas do 
meio e do fluido (YARIN et al., 2014). De acordo com o teorema, um grupo principal de 
variáveis deve ser considerado de forma a aparecer em todos os agrupamentos formados 
(WELTY, 2008). Tendo em vista a preocupação com as propriedades físico-químicas do 
fluido, as variáveis B”, y' e d' foram escolhidas para formar o grupo principal. Assim, o 
desenvolvimento dos grupos TT resulta em: 


quo spa adéixctsPemM sta Ul.ab 
L-ãa-b=-dr+D+i=0 (1.b) 
Masb=e+tdal]l=0 [1.6] 
Tid=b+ik-dl-1=0 Ud) 
ros Dd+ze+0=] =D (1. 


Resolvendo o sistema de equações composto pelas Equações 1.b-e de forma a 
compor um adimensional r,, tem-se que a=0; b=0,5; c=0,5; d=-1. Assim, o grupo 7, fica 


Dr E RA Ea ad a A RS aÃ ui Pa RD E 


, condado econsunis (1) 


famnad apare isa 


A obtenção dos demais grupos TT é análoga: 


Rea e pe ee me E e ÇÃO TUA (2.8) 


Le=Zg=b=304D+l= 0 (2.h) 
Matb=ec&d+D=0 (2.0) 
Tio=h4+2=72+0=0 (2d) 
[o 


0+0+dc+0+O Õ [2.6) 


densiiadestensão pap fesrincia agulha cplepor 


inter (21) 


e viscosidade 7 


Por possuírem a mesma dimensão que €, w'en' possuem desenvolvimento análogo. 
Assim: 
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Para os demais grupos: 


ao" peter ed gt =x a (5.a) 
Li =39=b=304043 0 (5.b) 
Masb=r+d+0=0 15.6) 
Tib=b+%-2- zh 5.) 
Avsosik+0+0=a (5.2) 
O na Mendo lento rpes itorato (51) 
ot appear egi= e Me x TExA! (6.a) 
L e=b=Ie4+0=]=D [:0+]] 
Masb-c+d+I=0 (8,6) 
TiDb=b4mk=7=1=0 [G.d] 
Meor0+2+0+0=0 [6.8] 
mf roca (8:) 
nin pêx exata e LEMA (7.a) 
L--3Ja-b-3e+0+1=0 (F.b) 
Mnth=cádraD=O [PC) 
ToD=b+4+Z-2- = (7d) 
Ao+D+ic+d+d=0 [7.e) 
n; = ftceeé deep ie pe [74] 


Como descrito por Welty (2008), os parâmetros adimensionais permitem relacionar 
a predominância das forças atuantes em uma superfície de controle durante um fenômeno 
físico. O parâmetro Tr, por exemplo, pode determinar a predominância entre as forças 
elétricas (campo elétrico, representado pela voltagem aplicada), viscosas (viscosidade) e 
de coesão entre as moléculas do fluido (tensão superficial) exercidas sobre o fluido. O 
parâmetro T1,, quando se considera a vazão de fluido na ponta de seção transversal circular 
da agulha, pode ser rearranjado como segue: 


sidade form Di agi EDITA O] ] rt eladade o rtri iate mulisar ira [iai ] 


= msi es : To ie ? | Piciialgadé mo erlris migra - u E Re [64] 


em que 0h é um adimensional conhecido como número de Ohnesorge, que relaciona 
exclusivamente as forças inerciais, viscosas e a tensão superficial do fluido (McKINLEY, 
2005) e Re é o número de Reynolds, bastante utilizado em Mecânica dos fluidos para 
a determinação de regimes de escoamento diversos e que relaciona forças inerciais e 
viscosas no fluido (WELTY, 2008). Neste caso em específico da Equação 8, no qual o 
comprimento característico foi tomado como sendo o diâmetro da agulha, a velocidade 


em questão é a razão entre a vazão do fluido polimérico e a área de seção transversal da 
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agulha. 
A combinação de adimensionais também pode resultar em outros parâmetros 
conhecidos. A razão Tr,/Tr, resulta em: 


E ELE E. — — [2] 


am | MEET 
idpe E mando duper E 


sendo a razão voltagem e distância entre agulha e coletor comumente designada 
de campo elétrico pseudo-homogêneo (RIEHLE, 1997; SHIN et al., 2001; GADKARI et al., 
2014; ZARGARIAN; HADDADI-ASL, 2017). No entanto, como a distância entre agulha e 
coletor está relacionada não apenas com a intensidade do campo elétrico, mas também 
com o tempo disponível para o desenvolvimento do jato e para a evaporação do solvente 
presente nas fibras durante a deposição (flight time) (TONG, WANG, 2010; ZAKARIA et al., 
2012; CAO et al., 2018), pode ser mais conveniente manter a distância agulha-coletor como 
variável separada da voltagem. 

Especialmente quanto à taxa de evaporação de solvente nas fibras durante sua 
deposição sobre o coletor, os adimensionais TT e Tr, juntamente a Tm, podem ser úteis 
para descrever o comportamento de secagem do solvente durante o processo, sendo o 
parâmetro 11, específico para coletores rotativos. Embora reconheça-se que as condições 
operacionais promovam diferentes taxas de evaporação e, com isto, a produção de fases 
cristalinas distintas no material formado (GILLER et al. 2010), trabalhos da literatura 
geralmente preocupam-se com o efeito da rotação do coletor sobre o alinhamento das 
fibras formadas, característica interessante para Engenharia de Tecidos (BOLAND et al., 
2001; MEDEIROS et al., 2008; TONG, WANG, 2010; ANINDYAJATI et al., 2015). Contudo, 
a literatura referente ao efeito da rotação de coletores rotativos sobre a fluidodinâmica do 
jato e sua secagem ainda é incipiente e estudos futuros podem ser realizados utilizando-se 
adimensionais como estes. 


31 METODOLOGIA 


Com o propósito de relacionar o conteúdo dos adimensionais com os resultados 
obtidos no presente trabalho, foram realizadas algumas predições utilizando o ferramental 
teórico para o cálculo da eficiência teórica de filtração utilizado por Oliveira et al. (2021). 
Também foram utilizados dados experimentais de Oliveira et al. (2021b) obtidos para 
as propriedades das soluções poliméricas de PVA, diâmetros de fibra e a densidade de 
empacotamento de meios filtrantes de nanofibras de PVA produzidas por electrospinning, 
exibidos na Tabela 2. Assim, se buscou estimar a eficiência de coleta, a queda de pressão 
e, consequentemente, o fator de qualidade avaliando faixas de tamanho de partículas 
(diferentes medianas) e tempos de fiação distintos dos utilizad s nestes trabalhos. 

Para a determinação da queda de pressão dos meios filtrantes limpos (ou queda de 


pressão inicial, que determina o consumo energético do processo), utilizou-se a equação 
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obtida por Davies (1953), que é expressa como: 


i 


gt Dr mo fl + nur" a] 


sendo H,ev respectivamente a viscosidade e a velocidade do ar. 

O fator de qualidade (g,) foi usado para avaliar a performance dos meios filtrantes 
baseado na eficiência de coleta global e na queda de pressão inicial, como segue (HINDS, 
1998): 


7 o Ut) 


em que E é a eficiência de coleta global calculada de acordo com a modelagem 
matemática de Oliveira et al. (2021a). 


C(% m/m) W (Pa.s) S (S/m) o (N/m) D (nm) a (=) 


0 3,27E+00 7,87E-02 4,89E-02 234 0,047 
0,25 4,91E+00 7,92E-02 3,65E-02 229 0,047 
0,50 3,54E+00 7,72E-02 3,43E-02 186 0,043 
0,75 3,64E+00 T,T6E-02 3,36E-02 191 0,052 
1,00 3,59E+00 7,66E-02 3,24E-02 180 0,054 


Legenda: C: concentração de surfactante; q: viscosidade; p: densidade; S: condutividade elétrica; o: 
tensão superficial; D: diâmetro de fibra; a: densidade de empacotament 


Tabela 2. Caracterização dos elementos filtrantes de Oliveira et al. (2021b) 


Como o estudo experimental de Oliveira et al. (2021b) consistiu no estudo da 
variação da concentração de surfactante Triton X-100 na solução polimérica sobre o 
diâmetro das nanofibras formadas, as propriedades dos meios filtrantes na presente análise 
poderiam ser ajustadas em função desta variável, a princípio. No entanto, isto limitaria a 
abrangência da análise, uma vez que os resultados obtidos não seriam representativos 
em estudos posteriores nos quais a concentração do polímero fosse a variável alterada 
na composição das soluções, por exemplo. Assim, se buscou relacionar as variáveis finai 
(eficiência, queda de pressão e fator de qualidade) em função das propriedades físicas 
das soluções. Para isto, o adimensional criado no estudo foi utilizado, uma vez que um 
único adimensional (17) reúne informações sobre a tensão superficial, a viscosidade e a 
condutividade elétrica. Dessa forma, a variável “concentração de surfactante” foi utilizada 
apenas para obter funções polinomiais simples que expressam a relação da concentração 
com as propriedades da solução a fim de se obter , como função destas. 

Salienta-se que os meios filtrantes produzidos experimentalmente por Oliveira et al. 
(2021b) foram produzidos sobre substratos constituídos de meios filtrantes de microfibra 
de celulose. No entanto, tais substratos foram utilizados apenas como suporte mecânico 


naquele trabalho, já que possuíam baixas eficiência e queda de pressão na filtraçã 
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de nanopartículas. Por isso, no presente trabalho desconsiderou-se a presença de tais 


substratos nos cálculos. 


41 RESULTADOS E DISCUSSÃO 


Para exemplificar a utilização dos adimensionais em dados reais, foram utilizados 
os dados da Tabela 2 para relacionar os adimensionais Tr, (referente às propriedades da 
solução polimérica) e 1, (contendo a variável diâmetro de fibra). Assim, a Figura 1 (a) 
apresenta os resultados calculados. A curva da Figura 1 (a) mostra um comportamento de 
máximo na concentração de 0,25% m/m de surfactante, o que está relacionado ao aumento 
expressivo da viscosidade em relação à solução sem surfactante (50% de variação desta 
propriedade frente às variações de 0,54% e -25% respectivamente para a condutividade 
elétrica e tensão superficial). Assim, a variação de Tr, para este caso está mais associada 
à variação da viscosidade e da tensão superficial, visto que o Triton X-100 utilizado é um 
surfactante não-iônico. O adimensional Tr, possui dependência inversamente proporcional 
com esta propriedade, conforme Equação 1.f. Em relação à dependência do diâmetro das 
fibras com as propriedades da solução, o comportamento exibido na Figura 1 (b) exibe 
decrescimento de acordo com a lei de potência (R? igual a 0,9708). 

Conforme mencionado anteriormente, a utilização da variável concentração de 
surfactante como variável independente limitaria a abrangência da análise teórica, uma vez 
que os resultados obtidos não poderiam não ser representativos em estudos posteriores 
nos quais a concentração do polímero fosse a variável alterada na composição das 
soluções, por exemplo. Dessa forma, a variável concentração de surfactante foi utilizada no 
presente estudo apenas para obter funções polinomiais simples que expressam a relação 
da concentração com as propriedades da solução a fim de se obter TT, como função destas, 
como mostra a Figura 2. Em outros estudos, outras variáveis (como a concentração de 
polímero) podem eventualmente ser utilizadas para obter os valores destas propriedades. 
Visto que a faixa de zero a 0,25% m/m apresentou valores de propriedades destoantes 
do restante da faixa de concentrações, possivelmente devido às diferentes interações do 
surfactante em solução (Oliveira et al., 2021b), optou-se por trabalhar apenas com a faixa 
de 0,25-1,0% m/m. A relação entre o adimensional Tr, e a concentração já está exibida na 


Figura 1 (a). 
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Figura 1. Adimensional TI, em função da concentração de PVA (a) e TI, em função de TT, (b) 


Com os ajustes polinomiais, foi possível obter pontos para serem plotados nos 
gráficos além daqueles quatro obtidos experimentalmente referentes às concentrações 
avaliadas (0,25; 0,50; 0,75; 1,0% m/m). Deste modo, para a obtenção das eficiência 
globais de coleta em função de Tm, foram utilizadas as equações do trabalho teórico 
de Oliveira et al. (2021a). A queda de pressão foi obtida a partir da equação de Davies 
(1953). Foram feitas simulações considerando medianas de distribuição em termos 
mássicos de 20; 50; 100 e 200 nm, com desvios padrões geométricos de 1,61 para todas 
as distribuições. Este último parâmetro é igual ao do trabalho experimental de Oliveira et 
al. (2021b), enquanto a mediana naquele trabalho era de -90 nm em termos mássicos. 
Além da faixa de concentração de surfactante (implícito no adimensional), foi considerado 
também o tempo de fiação. Manteve-se a velocidade de ar de 5 cm/s, bem como as demais 
condições ambientais e propriedades físicas da partícula, conforme utilizadas no trabalho 
de Oliveira et al. (2021b). 
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Figura 2. Propriedades das soluções de PVA utilizadas no trabalho de Oliveira et al. (2021b) em função 
da concentração de surfactante 


Para introduzir a discussão, a Figura 3 exibe como os dados se apresentam quando 
dados em termos do adimensional Tr, (gráficos à esquerda) e de concentração do surfactante 
(gráficos à direita), para a mediana de 20 nm 

Conforme se verifica na Figura 3, o comportamento das curvas com a variável 
concentração é bastante oscilatório. Os perfis utilizando o adimensional são mais suaves. 
Tal diferença de comportamentos decorre de o adimensional não ter uma relação linear com 
a concentração. Em outras palavras, o acréscimo entre 0,25 e 1,0% m/m de concentração 
não corresponde linearmente ao acréscimo entre 4,1 e -4,4x105 de TT. Isto ocorre uma vez 
que tal parâmetro aglutina várias propriedades da solução que possuem relações distintas 
com a concentração de surfactante (Oliveira et al., 2021b). No entanto, como exibido na 
Figura 3, os fenômenos físicos envolvidos na produção das fibras e sua relação intrínseca 
com as propriedades da solução continuam a ser descritos pelo adimensional. Além disso, 
a forma das curvas obtidas com Tr, permite uma parametrização mais simples (polinômios 
de baixo grau) para análises teóricas. Assim, os próximos gráficos (Figuras 4 e 5) serão 


dados apenas em termos do adimensional. 
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Figura 3. Eficiência (a, b), queda de pressão inicial (c, d) e fator de qualidade (e, f) em função de 


(gráficos à esquerda) e da concentração de surfactante (à direita 


Dos resultados das Figuras 4 e 5, fica claro que o fator de qualidade decai 


abruptamente com o aumento da mediana de tamanho de partícula utilizada. É preciso 


lembrar que o mínimo de eficiência fracionária (para cada diâmetro de partícula) de filtraçã 


encontra-se geralmente na região entre 100 e 200 nm (HINDS, 1998), o que prejudica a 


eficiência global de coleta. Visto que a queda de pressão é mais fortemente afetada pelo 


aumento da espessura da camada de fibras (que aumenta com o tempo de fiação) que 


a eficiência, o fator de qualidade decai bruscamente com essa variável, uma vez que a 


eficiência das fibras fiadas não está próxima o suficiente de 100% para que o fator de 
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qualidade aumente expressivamente. 

Uma consequência direta desta modelagem matemática é poder identifica 
previamente, a partir de uma queda de pressão limite, por exemplo, quais condições 
experimentais fornecem as melhores performances de filtração: obtendo-se as propriedades 
físicas do fluido, estando na faixa avaliada para o adimensional, seria possível calcular Tr, e 
realizar as devidas combinações com o tempo de fiação para obter a performance desejada. 
Estudos experimentais com as outras variáveis de operação do electrospinning (voltagem, 
distância agulha-coletor etc.) poderiam revelar as relações dos outros adimensionais com 
as variáveis de filtração e novas simulações poderiam ser feitas, assim por dia te. 

Aqui é preciso frisar que o ferramental teórico utilizado para as simulações (Oliveira 
et al., 2021a), superestimou os dados experimentais para partículas abaixo de 100 nm 
de diâmetro, seja pela simplicidade de algumas equações ou pela desconsideração de 
efeitos como a ressuspensão das partículas. Ainda assim, seria possível nortear pesquisas 
experimentais posteriores, que poderiam confirmar ou não o afastamento dos dados 
teóricos. 


51 CONCLUSÕES 


Neste estudo foram exploradas algumas aplicações dos adimensionais obtidos 
pelo teorema 11 de Vaschy-Buckingham para simplificar a análise da influência das 
propriedades físico-químicas sobre o processo de electrospinning. Comportamentos 
puderam ser evidenciados e sinalizaram que o diâmetro das fibras pode ser predito a partir 
das propriedades das soluções poliméricas (viscosidade, condutividade elétrica e tensão 
superficial), a depender de sua composição. Também foi possível utilizar o adimensional 
TT1 e os dados experimentais obtidos no presente trabalho para predizer a performance de 
meios filtrantes de PVA produzidos com diferentes tempos de fiação por electrospinning 


para diferentes faixas de tamanho de partícula. 
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Figura 4. Eficiência (a, d), queda de pressão inicial (b, e) e fator de qualidade (c, f) em função de , 
para mediana de 20 nm (a, b, c) e 50 nm (d, e, f) 
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Figura 5. Eficiência (a, d), queda de pressão inicial (b, e) e fator de qualidade (c, f) em função de , 
para mediana de 100 nm (a, b, c) e 200 nm (d, e, f) 
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O adimensional 171, atuando como intermediário entre as propriedades das soluções 
poliméricas e as variáveis de filtração, forneceu curvas de comportamentos mais suaves 
que utilizando-se a variável independente “concentração da solução”, podendo indicar 
maior praticidade e aplicabilidade para uso em análises teóricas deste tipo. 

Com a utilização dos adimensionais, futuros trabalhos poderão definir regimes de 
operação nos quais é possível obter fibras a partir de electrospinning de acordo com as 
propriedades da solução, bem como definir regimes de ocorrência de nódulos (beads) nas 
fibras de acordo com estas propriedades. Ainda, os resultados exibidos neste trabalho 
consideravam as condições operacionais constantes (voltagem, vazão de alimentação de 
solução polimérica, temperatura e umidade, rotação do coletor no caso de coletor rotativo 
etc). Entretanto, estes adimensionais podem ser utilizados em um planejamento estatístico 
de experimentos para a obtenção de correlações para o processo que independam do 
sistema de unidades adotado. 
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RESUMO: A economia de energia em processos 
industriais é um tema cada dia mais atual. 
Na indústria de alimentos, muita energia é 
empregada nos processos de esterilização de 
produtos alimentícios industrializados. Neste 
trabalho é apresentado o algoritmo de controle 
que foi desenvolvido e testado no aplicativo 
matlab/simulink"M visando permitir a economia 
de energia durante o processo de esterilização. O 
algoritmo consistiu em estimar a temperatura do 
produto alimentício utilizando modelo identificado 
para calcular o valor do índice de letalidade em 
tempo real. O modelo executa simultaneamente o 
cálculo da letalidade para a etapa de aquecimento 
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DE ALIMENTOS 


e para a etapa de resfriamento, a partir de um 
valor predeterminado da temperatura do produto. 
Estes valores são somados e o resultado é 
continuamente comparado com um valor de 
referência; quando a igualdade entre estes 
valores é verificada, um sinal discreto é gerado, 
permitindo interromper o aquecimento e iniciar 
o resfriamento. Os resultados obtidos foram 
validados por comparação com os disponíveis 
na literatura e revelaram boa aderência a estes 
últimos. 

PALAVRAS-CHAVE: Modelamento matemático, 
simulação computacional, esterilização industrial. 


INFERENTIAL CONTROL FOR FOOD 
STERILIZATION 
ABSTRACT: Energy saving in industrial 
processes is an increasingly current topic. In 
the food industry, a lot of energy is used in the 
sterilization processes of industrialized food 
products. This work presents the control algorithm 
that was developed and tested in the matlab/ 
simulinkTM application in order to allow energy 
savings during the sterilization process. The 
algorithm consisted of estimating the temperature 
of the food product using an identified model to 
calculate the value of the lethality index in real 
time. The model simultaneously performs the 
calculation of lethality for the heating step and 
for the cooling step, ffom a predetermined value 
of the product temperature. These values are 
summed and the result is continuously compared 
to a reference value; when the equality between 
these values is verified, a discrete signal is 
generated, allowing to stop heating and start 
cooling. The results obtained were validated by 
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comparison with those available in the literature and showed good adherence to the latter. 
KEYWORDS: Mathematical modeling, computer simulation, industrial sterilization. 


11 INTRODUÇÃO 


O processamento térmico de alimentos industrializados consiste no aquecimento 
de embalagens em autoclaves pressurizadas, que podem ser estacionárias ou rotativas, 
verticais ou horizontais, e/ou contínuas ou descontínuas. O processamento térmico deve 
ser especificado por especialistas que, para cada sistema formado por equipamento/ 
embalagem/alimento, definem as condições de processo, a saber: duração, temperatura 
de operação, temperatura mínima inicial do alimento, perfil de pressão e programa de 
desaeração da autoclave (Ohisson & Bengtsson, 2000). 

Em face destes requisitos de qualidade e de segurança, grande cuidado é tomado 
no cálculo destes tempos de processamento, assim como no controle do tempo e da 
temperatura vigentes durante os mesmos, para evitar que o alimento seja submetido a um 


subprocessamento ou a um sobreprocessamento. 


“Embora os processos possam ser criteriosamente definidos, a sua 
aplicação no dia-a-dia da produção industrial foge da idealidade. Problemas 
operacionais como a falta de energia, queda na geração de vapor de água 
e falha nos sistemas de controle ocorre com bastante frequência e causam 
desvios de processo pela flutuação da temperatura do meio de aquecimento. 
Outros desvios de processo podem ocorrer quando há alteração de 
parâmetros básicos de processo, tais como: temperatura inicial do produto, 
tempo de processo mais curto e temperatura do meio de aquecimento 
uniformemente mais baixa. Quando estes desvios acontecem, o lote - em geral 
- é reprocessado ou segregado para análise dos registros por um especialista 
que deve tomar decisões quanto à liberação deste ou o seu descarte. Estes 
procedimentos são dispendiosos, demorados e quando a opção é feita pelo 
reprocessamento, resultam num comprometimento da qualidade do alimento” 
(Germer, 1997). 


O objetivo deste trabalho é apresentar um modelo matemático inferencial que será 
empregado em conjunto com o sistema de controle da autoclave, com a finalidade de 
garantir determinado nível de letalidade previamente estabelecido para o microrganismo 
de referência. A letalidade é um índice tipicamente utilizado em processamento térmico 
de produtos alimentícios, usada para avaliar a mortalidade imposta ao microrganismo de 
referência. 

O emprego de modelos matemáticos capazes de predizer as características do 
produto em função de alterações nas condições de processo é fundamental para um bom 
desempenho dos sistemas de controle baseados em computador que sejam executados 


em tempo real (Akterian, 1999). 
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21 O PROCESSO DE ESTERILIZAÇÃO DE ALIMENTOS 


A esterilização térmica de alimentos enlatados em autoclaves pressurizadas é um 
dos métodos mais utilizados quando se diz respeito à preservação de alimentos. É um 
fator determinante ao bem-estar nutricional da população mundial. Eliminar esporos de 
bactérias, microrganismos nocivos ao corpo humano, é um dos principais objetivos na 
esterilização térmica. Esses microrganismos quando presentes nos alimentos produzem 
uma série de toxinas que mesmo sendo em quantidades infinitesim is, podem ser letais. 

Para isso existe o tratamento térmico desses produtos, que consiste em aquecer 
esses alimentos enlatados em autoclaves pressurizadas em temperatura especificada 
e por um período de tempo pré-definido, seguido por uma etapa de resfriamento. Este 
tratamento deve ser especificado de acordo com cada alimento e também depende muito 
do tipo de autoclave empregada. 

De um modo geral, a indústria de alimentos tem procurado utilizar de forma mais 
intensiva os sistemas de automação nos processos de enlatamento, com a finalidade de 
garantir que o processamento térmico imposto elimine, com segurança, os microrganismos 
indesejáveis no produto alimentício. Para se alcançar tal meta é essencial que o 
processamento térmico sofra um rígido controle toda vez que for empregado. Na Figura 
1 é apresentado um fluxograma de engenharia que mostra as malhas de controle de uma 
autoclave vertical típica. 

Para que a esterilização seja bem sucedida, é de fundamental importância o controle 
da temperatura à qual o alimento será exposto e do tempo de exposição àquela temperatura 
(Gedraite & Leonhardt, 1997). A malha de controle de temperatura da autoclave empregada 
neste trabalho é mostrada na Figura 1. Ela irá atuar após o período de exaustão, que tem 
por finalidade eliminar o ar presente no interior da autoclave. 

A cronometragem do período de exaustão e a do período de manutenção na 
temperatura necessária à letalidade dos microrganismos indesejáveis é feita por meio de 
um controlador lógico programável. 

Uma vez alcançada a esterilização desejada, deve-se proceder ao resfriamento 
da autoclave por meio da injeção de água. Nesta etapa do processo é muito importante 
o controle da pressão no interior da autoclave, para se evitar que os recipientes que 
contenham o alimento sejam danificados. A malha de controle de pressão, mostrada na 
Figura 1, irá atuar após o término da etapa de aquecimento, ajustando a vazão de ar 


comprimido que é alimentado à autoclave. 
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Figura 1 - Fluxograma de engenharia de uma autoclave estacionária vertical 


O controle da temperatura e da pressão no interior da autoclave é feito por 
controladores e indicadores eletrônicos microprocessados e o intertravamento da mesma 
é tipicamente realizado por meio de controlador lógico programável. Ele tem por finalidad 
principal impedir a colocação da mesma em operação se a tampa não estiver travada, além 
de impedir a abertura da tampa se a pressão interna for superior à pressão atmosférica. 
Para a medição da pressão e das temperaturas foram usados transmissores eletrônicos 
microprocessados, que permitem a obtenção de maior confiabilidade e maior precisão nas 
medições efetuadas (Gedraite & Leonhardt, 1997). 

Tipicamente, o sistema de controle da autoclave é empregado para gerar o perfil de 
temperatura desejado para a autoclave em função do tempo. O perfil de temperatura típico 
para uma autoclave estacionária é mostrado na Figura 2. Com base nesta figura pode-se 
distinguir, fundamentalmente, três fases de operação da autoclave, a saber: aquecimento, 
esterilização e resfriamento. A fase de aquecimento compreende as etapas de exaustão e 
de aquecimento propriamente dito da autoclave. A etapa de exaustão tem por finalidade a 
eliminação do ar existente no interior da autoclave, enquanto que a etapa de aquecimento 
visa atingir a temperatura pré-fixada para a esterilização. 

A primeira atividade a ser executada quando do início da operação da autoclave é a 
etapa de exaustão (ou remoção do ar). Considerando que o vapor de água é utilizado como 
meio de transferência de calor em autoclaves, devido, entre outras razões, aos maiores 
coeficientes de transferência de calor verificados, para que possamos tirar vantagem de 
seu uso, devemos eliminar o ar presente no interior da autoclave. A presença do ar, não 
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somente reduz a transferência de calor, mas, quando a exaustão não for total, o ar restante 
ao redor dos recipientes poderá agir como isolante e provocar a corrosão e queima de sua 
superfície (Rodrigues et al., 1998). 

Mesmo quando a autoclave se encontre com sua capacidade máxima ocupada, a 
maior parte do espaço total disponível é preenchida pelo ar. Em autoclaves horizontais, 
após seu preenchimento, 80% do seu espaço se encontra ocupado pelo ar. Em autoclaves 
verticais 60% da sua capacidade é tomada pelo ar. Para removermos todo este ar é 
necessário algum tempo de circulação do vapor com as válvulas de exaustão abertas, 
pois o ar deve ser desalojado dos espaços formados entre os recipientes. Por outro lado, 
devemos lembrar que, devido à superfície dos recipientes se encontrar em temperaturas 
inferiores à de saturação do vapor, ocorrerá uma condensação inicial. Por esta razão, 
a velocidade de fornecimento de vapor deve ser elevada para compensar as possíveis 
condensações, eliminar o ar presente e para que a temperatura de esterilização seja 
alcançada rapidamente (Gedraite, 1999). 

Uma vez que a temperatura recomendada para a exaustão seja atingida, deve- 
se manter a autoclave nesta temperatura pelo tempo necessário para se assegurar a 
eliminação do ar presente e, em seguida, providenciar que a válvula de exaustão seja 
fechada. Geralmente este intervalo de tempo é determinado experimentalmente e está 
baseado na concordância entre o valor da pressão lida no manômetro e a temperatura 
indicada no termômetro de mercúrio, ambos instalados no corpo da autoclave. Este 
procedimento não é uma indicação segura de que todo ar foi eliminado, no entanto, 
enquanto esta concordância não for verificada, deve-se prosseguir com a desaeração 
(Leonhardt, 2003; Gedraite, 1999). 
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Figura 2 - Perfil de temperatura da autoclave em função do tempo 


(adaptado de Leonhardt, 2003) 
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A seguir, deve-se continuar o fornecimento de vapor de água, até que seja 
alcançada a temperatura adequada para a esterilização do produto alimentício. A partir 
deste ponto, deve-se iniciar o fechamento gradual da válvula de controle que é empregada 
para ajustar o valor da vazão de vapor de aquecimento. Este procedimento evita uma 
diminuição acentuada da temperatura, que geralmente ocorre quando esta válvula é 
fechada abruptamente. A partir deste ponto, tem início a cronometragem do tempo de 
processamento, o que é feito pelo sistema de automação da autoclave. Durante a fase de 


esterilização, os seguintes cuidados são necessários (Toledo, 1994): 


* Manter a temperatura da autoclave um grau acima da temperatura recomenda- 
da, para compensar flutuações 


* Verificar periodicamente, através do termômetro de mercúrio, a temperatura da 
autoclave; 


* Verificar se, em todos os sangradores, o vapor está escoando normalmente 


* verificar o tempo de processamento e, uma vez esgotado este tempo, iniciar 
imediatamente a operação de resfriamento. 


O resfriamento do sistema é conduzido na própria autoclave, sob pressão controlada. 
Esta pressão pode ser mantida praticamente constante, tanto com a introdução de vapor 
de água quanto com a introdução de ar comprimido. Quando utilizamos vapor de água para 
obtermos a pressão desejada, deve-se adotar os seguintes procedimentos (Leonhardt, 
2003; Gedraite, 1999): 


1. Decorrido o tempo pré-fixado para a esterilização, deve-se providenciar a 
interrupção do fornecimento de vapor de água e o fechamento da válvula de dreno. 
A elevação da pressão da autoclave é feita pela alimentação de vapor de água, 
fornecido pelo topo, de maneira a se manter uma pressão cerca de 2 psi maior que 
a pressão de esterilização. 


2. Atingida a pressão desejada, introduzimos água e vapor pela base até termos 
preenchido um quarto de capacidade da autoclave com água em ebulição. Isto feito, 
fechamos o vapor da base e, controlando a vazão de vapor no topo, continuamos 
a introduzir água. Este procedimento é necessário para evitarmos a condensação 
súbita do vapor, o que provocaria uma deformação permanente dos recipientes. 


3. Quando a água chega ao topo, o sistema de automação faz com que a válvula de 
extravasamento seja gradualmente aberta, executando assim o controle da pressão, 
o que deve ser feito até que a temperatura do produto atinja 60 -70 ºC. Devemos ter 
em mente que, se o escoamento através da canalização de extravasamento for bem 
controlado, a pressão de resfriamento pode ser mantida sem a introdução de vapor 
de água no equipamento. 


4. Tendo alcançado esta temperatura, deve-se iniciar o relaxamento da pressão, que 
dependerá, fundamentalmente, do tipo de produto que está sendo processado. Se 
o alimento é aquecido por convecção, o relaxamento deve ser feito na base de 2 
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psi/minuto. Caso seja aquecido por condução, a redução de pressão deve ser mais 
lenta e a velocidade de relaxação deve ser de 1 psi/minuto. Dependendo do tipo de 
produto, pode ser necessário manter a pressão de resfriamento por 10 a 15 minutos 
após a temperatura do produto ter atingido 60 - 70ºC. 


5. Uma vez atingida pressão manométrica nula, o alimento já se encontra praticamente 
frio. Porém, caso seja necessário, devemos prosseguir com o resfriamento até que a 
temperatura do produto atinja 35 - 40 ºC. Esta temperatura final é conveniente para 
favorecer a evaporação da água superficial dos recipientes, quando os retiramos da 
autoclave. 

Quando a pressão na operação é conseguida através de ar comprimido, o 
procedimento a ser adotado é praticamente o mesmo. Terminada a esterilização, interrompe- 
se o fornecimento do vapor de água e fecha-se a válvula de dreno. Tem início a injeção de 
ar comprimido pelo topo da autoclave, de modo a se obter a pressão desejada, após o que 
tem início o fornecimento de água de resfriamento pela base. Naturalmente, a tendência 
da pressão é cair, daí a necessidade de um controle adequado da pressão interna da 
autoclave. No instante em que a água atinge o topo, o sistema de automação abre a válvula 
de dreno, mantendo a pressão estipulada para o resfriamento. O relaxamento da pressão, 
que será realizado de acordo com o tipo de produto, só deve se iniciar após a temperatura 
do produto ter atingido 60 - 70ºC. 

Tanto no caso da pressão de resfriamento ser obtida com vapor, como com ar 
comprimido, a vazão de água, bem como sua circulação, deve satisfazer a certos requisitos 
básicos como: (i)- para encher a autoclave devemos gastar, no máximo, 5 minutos e 
(ii)- no caso da autoclave ter entrada de água pela base e pelo topo, deve-se inverter 
a introdução de água, ou seja, inicialmente a água é fornecida durante 5 minutos pela 
base, em seguida durante mais 5 minutos pelo topo, e assim alternadamente até o final do 
processo (Leonhardt, 2003). 

A preocupação em se evitar o desperdício do agente de aquecimento utilizado, em 
geral o vapor de água saturado, através da otimização do tempo de processamento térmico 
imposto ao produto alimentício, assegurando que o nível de letalidade acumulada desejado 
seja alcançado, motivaram o estudo da possibilidade de aplicação de uma estratégia de 
controle mais elaborada, que será apresentada na sequência. 

O sistema de controle considerado neste trabalho foi baseado no emprego de 
modelo matemático a parâmetros concentrados para estimar a temperatura do produto 
alimentício em tempo real. O modelo utilizado foi do tipo semiemprírico identificado, obtido a 
partir de experimentos realizados com o sistema em estudo (Gontijo et al, 2021; Larousse 
& Brown, 1997). 
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81 CONSIDERAÇÕES ACERCA DA ESTRATÉGIA DE CONTROLE DO 
PROCESSO 

O sistema de controle da autoclave é empregado visando gerar o perfil de temperatura 
desejado para a autoclave, conforme já abordado em Veloso et al. (2015). 

A esterilização de um produto alimentício depende essencialmente da taxa de 
penetração de calor no mesmo (Teixeira & Balaban, 1996). 

Quando produtos alimentícios são colocados numa embalagem, que por sua vez é 
colocada numa autoclave e processada termicamente através da condensação do vapor 
de água saturado, a taxa de transferência de calor para os mesmos depende dos seguintes 
fatores (Fontana et al., 2001): 

* | coeficientes de transferência de calor 

* propriedades físicas do produto e da embalagem; 

* | diferença de temperatura entre o produto alimentício e autoclave; 
* | tamanho da embalagem. 

Neste trabalho optou-se por utilizar um produto sólido cujo mecanismo de transporte 
de calor seja essencialmente condutivo. 

A principal vantagem advinda do uso de modelos matemáticos capazes de estimar o 
histórico da temperatura do produto é o fato de que se pode analisar flutuações imprevistas 
na temperatura da autoclave. Assim sendo, na ocorrência de uma diminuição inesperada 
da temperatura da autoclave, o modelo fornecerá rapidamente o histórico das temperaturas 
e calculará o valor do nível de letalidade atingido. 

Assim, na ocorrência de uma perturbação na temperatura da autoclave, o modelo 
deverá prever a resposta correta da temperatura no interior do produto para este desvio 
(Larousse & Brown, 1997). 

Para que haja consistência entre o modelo apresentado neste trabalho e a realidade, 
algumas considerações foram feitas (Gontijo et al., 2021). As mesmas são apresentadas 
a seguir: 

. modelo à parâmetros concentrados; 
*. o produto em estado sólido; 


* | otransporte de calor no interior do produto ocorre por condução; 


* | o transporte de calor na interface superfície da embalagem/vapor de água de 
aquecimento ocorre por convecção; 


* | o headspace foi considerado desprezível no interior da embalagem. 


O modelo matemático foi desenvolvido tomando-se por base a aplicação de um 
balanço global de energia no volume de controle correspondente ao conjunto embalagem/ 
produto (Teixeira et al., 1999). A equação do balanço global de energia aplicada ao sistema 
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em estudo é apresentada na Eq. (1). 
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O calor trocado por unidade de tempo entre o fluido de aquecimento/resfriamento 
utilizado no interior da autoclave e o conjunto embalagem/produto alimentício foi modelado 
com base na Eq. (2). 


O = UA[TÇE) = Tilt] (2) 


Considerando a representação matemática do comportamento da temperatura Ta 
em função do tempo, na etapa de aquecimento, pode-se escrever a Eq. (3). 
O m HATE) = To(t)] (2] 
Isolando-se os termos que contém a derivada na Eg. (3), resulta a Eq. (4). 
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A simulação foi conduzida utilizando-se o modelo identificado do processo, o qual 


foi baseado na Eq. (4). Na Figura 3 é apresentado o diagrama de simulação desenvolvido 
neste trabalho. 
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Figura 3 - Diagrama de simulação do processo e sistema de controle 


O modelo apresentado na Figura 3 pode ser subdividido em três partes distintas, 
a saber: (i)- o modelo que representa a etapa de aquecimento, marcado na figura com o 
número 1, (ii)- o modelo que representa a etapa de resfriamento, marcado na figura com 
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número 2 e (iii)- o modelo que representa a estratégia de controle inferencial propriamente 
dita, marcado na figura com o número 3 

O modelo que representa a etapa de aquecimento foi explicado em Veloso et al. 
(2015), tendo sido no presente modelo incluída a equação utilizada para o cálculo do índice 
de letalidade (F,) correspondente à etapa de aquecimento. 

A partir do comportamento simulado da temperatura da autoclave, considerando uma 
perturbação do tipo degrau em seu valor, é calculada a temperatura do produto alimentício 
em estudo, utilizando o conjunto de funções (blocos) identificados por MODELO Talim. 
Com base no valor da temperatura do alimento é calculado o valor de F, correspondente 
à etapa de aquecimento, por meio do bloco de função identificado por Calc F, aquecim. O 
valor de F, calculado é enviado para os blocos que representam a estratégia de controle 
inferencial do processo. 

Quando a temperatura do produto alimentício se torna maior do que 100ºC, a 
estratégia de controle inferencial é acionada. Isto ocorre quando a temperatura do produto 
alimentício ultrapassa em 20ºC a temperatura anteriormente citada. O critério utilizado para 
a escolha do valor de 100ºC baseia-se no fato de que o valor de F, não é significativo para 
valores menores de temperatura (Gedraite et al., 1998). 

A partir da condição citada no parágrafo anterior, é iniciada a simulação da etapa de 
resfriamento, representada no modelo por meio do chaveamento do valor nulo para o valor 
120ºC, que é a diferença de temperatura aplicada no degrau da etapa de aquecimento. O 
conjunto de blocos identificados por MODELO Talim calculam a temperatura do produto 
alimentício e o bloco de função identificado por Calc F, resfriam calcula o valor de F, 
correspondente à etapa de resfriamento. 

Os valores calculados de F, nas etapas de aquecimento e de resfriamento são 
somados, dando origem ao valor do F, total, o qual será comparado com o valor de F, target 
para calcular o sinal de desvio que será utilizado pela estratégia de controle inferencial para 
interromper o processo. Este será interrompido quando o valor do sinal de desvio for menor 
ou igual a um limite preestabelecido pelo operador do processo. 

Avalidação do modelo utilizado para o cálculo da temperatura do produto alimentício 
foi feita pela comparação dos resultados obtidos pela simulação computacional com 
aqueles apresentados na literatura (Gedraite et al., 1998). 


41 RESULTADOS & DISCUSSÃO 


Quando o sistema de controle inferencial é posto a funcionar, ele continuamente 
compara o valor de F, target com o valor de F, total para interromper o aquecimento e 
iniciar o resfriamento. O cálculo do valor de F, total é feito com base no valor simulado de 
F, resfriam, o qual começa a ser computado quando a temperatura do alimento ultrapassa 


em 20ºC o valor estabelecido como critério que é 100ºC. Na Figura 4 é apresentado o 
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comportamento temporal da temperatura do alimento e do sinal de controle considerado, 
na qual se pode perceber que aproximadamente no instante t = 70 minutos a temperatura 
do produto alimentício atende ao critério estabelecido. 

O sistema de controle inferencial passa então a computar a contribuição do 
resfriamento no cálculo do valor de F, total. O sistema de controle continua atuando até 
que a diferença entre os valores de F, total e de F, target seja nula. Quando isto acontece, 
aparece na tela do computador a mensagem gerada pelo bloco de função Check Static 
Gap informando que a condição definida para interromper o processo de aquecimento foi 


alcançada. Esta condição é ilustrada na Figura 5. 
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Figura 4 - Condição de início do cálculo do valor de F, resfriam 


No processo estudado neste trabalho, o instante em que o valor de F, total iguala o 
valor de F, target é igual a 202,98 segundos. Quando esta condição é atingida, a estratégia 
de controle inferencial interrompe a execução da etapa de aquecimento e inicia a etapa de 
resfriamento, como mostrado na Figura 6. O bloco de função do aplicativo simulink"M que 
é responsável por esta ação é o bloco seletor lógico configurado no início de resfriamento. 
O papel deste bloco é comparar o sinal gerado pela diferença entre o valor de F.(t) 
integrado ao longo do tempo e o valor previamente estabelecido de 65 com o valor definid 
internamente àquele bloco como ponto de ajuste. Quando esta condição é atendida, o bloco 
seletor lógico início de resfriamento passa a gerar o sinal correspondente à diminuição 
da temperatura de 120ºC para 20ºC, iniciando a simulação da etapa de resfriamento. Na 
Figura 7 é apresentado o comportamento temporal do valor calculado de F, aquecim e na 


Figura 8 o comportamento temporal do valor calculado de F, resfriam. 
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Figura 5 - mensagem informando que a condição definida para a interrupção do processo foi alcançad 
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Figura 6 - Condição de início do cálculo do valor de F, resfriam e de efetivo início da etapa de 
resfriamento 
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Figura 7 - Comportamento temporal do valor calculado de F, aquecim 
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Figura 8 - Comportamento temporal do valor calculado de F, resfriam 


51 CONCLUSÕES 


A estratégia de controle inferencial proposta neste trabalho foi considerada 
satisfatória para fins de controle do processo de esterilização de produtos alimentícios 
empregando o calor como agente de esterilização. 

Pode-se, portanto, empregar o modelo matemático em estudo para fins de controle 
regulatório de processo, visando melhorar o desempenho da operação da autoclave. 

Por se tratar de modelo a parâmetros concentrados, o modelo matemático 
desenvolvido poderá ser facilmente implementado em código computacional que permita 
a troca de informações em tempo real com o sistema de controle instalado em autoclaves 
estacionárias tipicamente utilizadas na indústria alimentícia. 

Este modelo pode ser adaptado para estudo do transporte de calor em diferentes 
embalagens. 


61 NOMENCLATURA 


A= Área de troca térmica da embalagem (m?) 
Copa Calor específico médio do produto alimentício (J/ kg - ºC) 
Coeficiente global de transporte de calor entre o fluido e a embalagem d 


Ee produto alimentício (W/ m? - ºC) 

= Massa do produto alimentício (kg) 

= Calor transferido por unidade de tempo (W) 
T(D= Temperatura do produto alimentício no instante considerado (ºC) 
T(B= Temperatura da autoclave no instante considerado (ºC) 


= Variação da temperatura do produto alimentício com o tempo (ºC/s) 
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RESUMO: Com uma extensa gama de aplicações 
industriais, o ácido cítrico desperta o interesse 
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DESCONTÍNUO 


no estudo de sua biossíntese e parâmetros 
processuais produtivos. Atualmente, (o) 
Aspergillus niger é o principal fungo para o cultivo 
cítrico. Neste cenário, o presente estudo avalia a 
produção cítrica, em sistema descontínuo, pelo 
Aspergillus awamori, investigando a influência da 
temperatura no meio. Assim, foram conduzidos 
processos em três temperaturas à 30, 33 e 35ºC 
apresentando uma concentração final cítrica 
de 9, 6 e 3 g/L, respectivamente. Relativo ao 
consumo de substrato, para os processos a 30, 
33 e 35ºC obteve-se reduções de 45, 45 e 85 
g/L da alimentação original, nesta ordem. Em se 
tratando do crescimento fúngico, os processos a 
30, 33 e 35ºC apresentaram 37, 32 e 30 g/L em 
sua concentração máxima. Além disso, avaliou- 
se o pH e a acidez total no meio de cultivo. 
PALAVRAS-CHAVE: Ácido Cítrico, Aspergillus 
awamori, Aspergillus niger, bioprocessos, 
fermentação. 


PRODUCTION OF CITRIC ACID 
BY ASPERGILLUS AWAMOR!I AND 
INVESTIGATION OF TEMPERATURE 
EFFECTS IN DISCONTINUOUS SYSTEM 
ABSTRACT: With an extensive range of industrial 
applications, citric acid arouses interest in the 
study of its biosynthesis and production process 
parameters. Currently, Aspergillus niger is the 
main fungus for citrus cultivation. In this scenario, 
the present study evaluates citrus production, in 
a discontinuous system, by Aspergillus awamori, 
investigating the influence of temperature on the 
medium. Thus, processes were carried out at 
three different temperatures, 30, 33 and 35ºC, 
presenting a final citrus concentration of 9, 6 and 
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3 g/L, respectively. Regarding the consumption of substrate, for the processes at 30, 33 and 
35ºC, reductions of 45, 45 and 85 g/L of the original feed were obtained, respectively. In 
terms of fungal growth, the processes at 30, 33 and 35ºC presented 37, 32 and 30 g/L at 
their maximum concentration. In addition, pH and total acidity in the culture medium were 
evaluated. 

KEYWORDS: Citric Acid, Aspergillus awamori, Aspergillus niger, bioprocesses, fermentation. 


11 INTRODUÇÃO 


Considerado um ácido orgânico essencial e multifuncional, o 2-hidroxopropano- 
1,2,3-tricarboxílico (Figura 1), de nome usual ácido cítrico, possui tanto uso em aplicações 
domésticas, bem como industriais, haja vista ser uma molécula versátil, a qual apresenta 
propriedades acidulantes, flavorizantes, antioxidantes, quelantes e plasticizantes, conforme 
Nigan (2009). 


e. oo A “0H 


Figura 1: Fórmula estrutural do ácido cítrico. 


Fonte: Barros (2018) 


Primeiramente, cabe referir que apesar do processo de extração do ácido cítrico, 
por meio de frutas cítricas, haver dominado o esquema produtivo durante o século XX, na 
indústria contemporânea há uma ênfase no sistema de produção de cultura submersa, o 
qual dispõe de processos biotecnológicos, com utilização de microrganismos, normalmente 
do gênero Aspergillus (Figura 2) e que segundo Soccol et al. (2006), apresenta rendimento 


significativamente superior a aqueles obtidos através de síntes química. 


fai (bi 


Figura 2: (a) conidióforo de Aspergillus visualizado através de microscopia. (b) Esquema representativo 
da estrutura reprodutiva do fungo Aspergillus. 
Fonte: Adaptado de Abreu (2012). 
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Relativo à excreção cítrica, este fenômeno ocorre em inúmeros microrganismos, 
inclusive em bactérias, fungos e leveduras. Classificado como um composto intermediário 
do metabolismo normal, seu acúmulo é condicionado a um desequilíbrio metabólico, o 
qual pode ser induzido por fatores externos (temperatura, agitação, etc...), no entanto, de 
acordo com Kubicek e Rôhr (1986), poucas espécies são capazes de secretar quantidades 
apreciáveis do ácido cítrico. 

Neste âmbito, a avaliação de novas espécies para produção de ácido cítrico é 
de grande relevância industrial, devido à crescente demanda e vasta aplicabilidade. 
Outrossim, segundo Goldberg et al. (2006) e Rymowicz et al. (2008) o volume desta 
produção é crescente, a uma taxa que excede a de qualquer outro ácido orgânico obtido 
por microrganismo. Destarte, o presente artigo expõe um estudo da utilização do fungo 
Aspergillus awamori, em condições variadas de temperatura, com enfoque na produção de 
ácido cítrico. Ademais, traz-se uma avaliação sumária da relação entre o parâmetro físico 


investigado com o crescimento celular, consumo de substrato e o pH. 


21 RESULTADOS E DISCUSSÃO 


2.1 Efeitos da Temperatura na Produção Cítrica por Aspergillus awamori 


O microrganismo utilizado foi o A. awamori NRAL 3112, obtido primariamente na forma 
liofilizada, e posteriormente submetido à preparação do inóculo. O processo fermentativo 
foi conduzido no biorreator New Brunswick Scientific modelo Bio Flo Ill. Afim de precisar 
os açúcares redutores totais, fora empregado o método DNS modificado, desenvolvido por 
Zanin e Moraes (1987). Para aferir a concentração celular utilizou-se o método descrito por 
Olivo (1998). A determinação do ácido cítrico foi realizada em cromatógrafo de fase líquida 
de alta eficiência (High Pressure Liquid Chromatography, HPLC, Varian 920-LC), sendo 
operado a coluna analítica Aminex HPX-87Hº (300 mm x 7,8 mm) fabricada pela BIORADS. 
A acidez total foi analisada segundo o método descrito pelo Instituto Adolfo Lutz (2008). 
Para mais, detalhes da metodologia de preparo do inóculo, meio de cultivo, biorreator, 
materiais e métodos específicos podem ser encontrados no trabal o de Barros (2018). 

Deste modo, com o propósito de se avaliar a influência da temperatura no processo 
produtivo de ácido cítrico, realizaram-se processos nas temperaturas de 30, 33 e 35 
ºC (Figura 3), sob uma condição inicial de substrato de 200 g/L, agitação de 200 rpm e 
inóculo em 10% do volume inicial do meio a ser cultivado. O tempo total de cultivo foi de, 


aproximadamente, 300 horas. De maneira que, variou-se apenas a temperatura. 
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Figura 3 — (A) Produção de ácido cítrico ([ÁC] (m)) por A. awamori, consumo de glicose (S (+), 
produção de biomassa fúngica (X (4)) e (B) correlação entre o ácido cítrico (4), ácidos totais 
(») e os valores de pH (nm). Sendo que, os processos 1, 2 e 3 foram conduzidos à 30, 33 e 35ºC, 
respectivamente. 


Referente ao crescimento fúngico, observa-se na Figura 3 (A) que a partir das 168 
horas de processo, há uma queda na curva de crescimento decorrente da diminuição 
da multiplicação celular, embora o biorreator permanecesse em condição de agitação 
constante, ocorreu acúmulo de células nas zonas mortas e superfície do equipamento, de 
forma a “reduzir” a concentração celular dispersa no meio. Assim a concentração celular no 
ponto de máxima desta, isto é, aproximadamente às 168 horas, a 30ºC foi de 37 g/L, Figura 
3 (A1),a 33ºC de 32 g/L, Figura 3 (AZ) e a 35ºC de 30 g/L, Figura 3 (A3). O comportamento 
do A. awamori mostrou que, próximo a 30ºC, Figura 3 (A1), houve um crescimento 
celular mais acentuado do que em 33 e 35ºC, Figura 3 (A2) e (A3), contrariando o que foi 
identificado por Kolicheski (1995), o qual em seus estudos afirmou que, o cultivo de células 
fúngicas, em específico do A. niger, para a formação do ácido cítrico em temperaturas mais 
elevadas ocorria rapidamente e com abundante crescimento de micélios. 

Relativo ao consumo de substrato, os processos à 30 e 33ºC, Figura 3 (A1) e 
Figura 3 (A2) respectivamente, apresentaram semelhanças nas tendências de suas curvas 
e redução em suas fontes carbonáceas de 45 g/L, ou seja, 22% do montante original. 
Diferentemente, o processo à 35ºC, Figura 3 (A3), mostrou maior consumo de substrato, 
da ordem de 85 g/L, isto é, 43% do total alimentado. Desta forma, o 4. awamorie o A. 
niger, em temperaturas mais elevadas denotaram uma similaridade comportamental, 
face ao estudado por Kolicheski (1995), o qual refere que a produção cítrica, por 4. 


niger, em temperaturas mais elevadas provoca uma alta taxa de oxidação do açúcar e 
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consequentemente baixa produção do ácido cítrico. 

Com relação à produção cítrica, verifica-se que ambas as curvas delineiam um perfi 
na forma “S” sigmoidal crescente repetida, Figura 3 (A), com dois intervalos, os quais são: 
primeiro de adaptação, com baixa geração de produtos e o segundo, no qual percebeu-se 
um incremento significativo na produtividade do ácido cítrico e consequente redução no 
pH do meio, chegando a cerca de 3, valor considerado adequado por Pandey et al. (2016). 
Desta maneira, nos processos à 30ºC a concentração cítrica final foi de 9 g/L, Figura 3 (A1), 
à 33ºC uma produção de 6 g/L Figura 3 (A2) e a 35ºC obteve-se 3 g/L Figura 3 (A3), sendo 
que a relação percentual da produção entre o processo à 30ºC e à 35ºC é de 77% a mais 
para o primeiro. 

Correspondente aos ácidos totais, em todos os processos, mostrados na Figura 
3 (B), suas curvas geraram perfis que se assemelham aos da produção do ácido cítrico. 
Contudo, destaca-se o processo à 35ºC Figura 3 (B3), no qual houve maior discrepância 
entre a curva de ácido cítrico e ácidos totais, o que conforme relatado por Sodeck et al. 
(1981) ocorreu pois, em temperaturas que excedem o limite máximo, observa-se à geração 


de outros ácidos do ciclo do ácido tricarboxílico, como o oxálico, na produção cítrica. 


31 CONCLUSÃO 


A respeito do crescimento fúngico, o aumento da temperatura, acima de 30ºC, 
desfavorece o desenvolvimento celular no meio de cultivo. Atinente ao consumo de 
substrato, ocorre significativo aumento deste com a elevação da temperatura, dentro do 
espectro estudado, e em consequência há o surgimento de outros ácidos do ciclo do ácido 
tricarboxílico. Em consideração a produção cítrica, pode-se inferir que, em temperaturas mais 
amenas, próximas de 30ºC, há melhores condições de formação do produto. Depreende- 
se que, em análise comparativa entre 4. awamori e o A. niger, há necessidade de estudo 
mais abrangente, para que o microrganismo em questão seja competitivo industrialmente. 
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RESUMO: O Eucalyptus é um dos gêneros 
florestais mais cultivados no Brasil e sua 
produção origina diversos resíduos como raízes, 
galhos e folhas. Uma forma de reaproveitamento 
desse material é aplicando pré-tratamentos 
químicos e biológicos para fracionar seus 
componentes lignocelulósicos, e assim ser 
destinados para a produção de diversos 
subprodutos. O objetivo deste trabalho foi 
estudar diferentes pré-tratamentos para madeira 
in natura e madeira tratada termicamente a 150º 
C do gênero Eucalyptus, a fim de atingir uma 
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maior deslignifi ação e preservar o máximo de 
carboidratos. No estudo deste trabalho foram 
realizados testes em que o Eucalyptus in natura 
foi submetido a hidrólise com ácido peracético 
(PAA) em temperaturas de 60ºC e 90ºC, 
diferindo o tempo de contato. Sendo assim, 
observou-se que a elevação da temperatura 
atuou como um dos fatores fundamentais, visto 
que promoveu uma maior remoção da lignina 
(43%). Consequentemente, essa deslignificação 
implicou positivamente em maior extração de 
xilose (40%) e glicose (15%) com maior carga 
de enzimas 22,5 FPU (Unidades de Papel 
Filtro). Outro pré-tratamento foi realizado com 
eucalyptus tratado termicamente a 150º C, na 
presença de hidróxido de sódio e posteriormente 
tratado com enzima lacase para remoção da 
lignina e frações da hemicelulose. Observou-se 
que a amostra contendo lacase com adição de 
substrato ABTS apresentou maior intensidade 
na espectrofotometria UV/VIS e com isso maior 
extração de xilose (37%) e glicose (10%). Diante 
destes dados, o pré-tratamento com ácido 
peracético a 90 “C apresentou maior extração de 
xilose e glicose, quando comparado com o pré- 
tratamento com NaOH e lacase. 
PALAVRAS-CHAVE:  Hidrólise enzimática. 
Remoção de lignina. Pré-tratamento. Extração 
de açúcares. 


EXTRACTION OF XYLOSE AND 
GLUCOSE FROM EUCLYPTUS sp. PRE- 
TREATED WITH PERACETIC ACID AND 

SODIUM HYDROXIDE 
ABSTRACT: Eucalyptus is one of the most 
cultivated forest genera in Brazil and its 
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production originates various residues such as roots, branches and leaves. One way to 
reuse this material is by applying chemical and biological pre-treatments to fractionate its 
lignocellulosic components, and thus be destined for the production of various by-products. 
The objective of this work was to study different pre-treatments for wood in natura and wood 
heat treated at 150º C of the genus Eucalyptus, in order to achieve greater delignification 
and preserve the maximum amount of carbohydrates. In the study of this work, tests were 
carried out in which Eucalyptus in natura was subjected to hydrolysis with peracetic acid 
(PAA) at temperatures of 60ºC and 90ºC, differing the contact time. Thus, it was observed 
that the increase in temperature acted as one of the fundamental factors, since it promoted a 
greater removal of lignin (43%). Consequently, this delignification positively implied in greater 
extraction of xylose (40%) and glucose (15%) with a greater load of enzymes 22.5 FPU (Filter 
Paper Units). Another pre-treatment was carried out with eucalyptus heat treated at 150ºC in 
the presence of sodium hydroxide and subsequently treated with laccase enzyme to remove 
lignin and hemicellulose fractions. It was observed that the sample containing laccase with 
the addition of ABTS substrate showed greater intensity in the UV/VIS spectrophotometry 
and with that, greater extraction of xylose (37%) and glucose (10%). In view of these data, 
the pretreatment with peracetic acid at 90 “C showed higher extraction of xylose and glucose, 
when compared with the pretreatment with NaOH and laccase. 

KEYWORDS: Enzymatic hydrolysis. Lignin removal. Pre-treatment. Sugar extraction. 


11 INTRODUÇÃO 


O Brasil tem uma imensa diversidade de biomassas lignocelulósicas apresentando 
um potencial para aproveitamento energético e químico em larga escala, a exemplo dos 
resíduos gerados durante a exploração florestal e nas serralherias (MACÊDO, 2012). 
Biomassa celulósica, derivada de fontes que não são alimentos, tal como árvores, está 
sendo utilizada como matéria-prima para produção de etanol celulósico. Na parede celular 
vegetal, celulose é encontrada em combinação com lignina e hemicelulose resultando em 
um sólido complexo e recalcitrante chamado de lignocelulose (CASTOLDI et al., 2014). 

Esses materiais possuem uma rede complexa e resistente formada principalmente 
por lignina (10% a 30%), hemicelulose (15% a 35%) e celulose (30% a 50%) e as 
concentrações de cada um desses elementos variam de acordo com o tipo de matéria-prima, 
idade e estágio vegetativo (JEFFRIES, 1990). Os resíduos de materiais lignocelulósicos 
podem ser aplicados na produção de alimentos, combustíveis, insumos químicos, enzimas 
e bens de consumo diversos (HAMELINCK, 2005). 

Materiais baseados em açúcar e amido tais como cana-de-açúcar e grãos são dois 
grupos de matérias-primas atualmente usados como fontes principais para produção de 
etanol. O terceiro grupo é os materiais lignocelulósicos representando a opção mais viável 
para produção de etanol (TALEBNIA et al., 2010). No intuito de produzir etanol celulósico, 
diferentes estudos têm sido desenvolvidos, como sacarificação e fermentação simultâneos 
(SSF), sacarificação simultânea à co-fermentação (SSCF), hidrólise e fermentação 


simultâneos (SHF), hidrólise e co-fermentação das hexoses com pentoses (SHCF) e 
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bioprocessamento consolidado (CBP) (SILVEIRA; VANELLI; CHANDEL, 2018). 

A estrutura dos materiais lignocelulósicos os torna resistentes à degradação tanto 
dos polissacarídeos quanto da lignina (BRIENZO, 2010). Para o fracionamento da biomassa 
lignocelulósica é necessário a aplicação de tratamentos físicos, químicos e biológicos ou 
em conjunto. O pré-tratamento é uma forma que visa a diminuição da recalcitrância dos 
materiais para facilitar o fracionamento. Tratamentos químicos em condições alcalinas e 
altas temperaturas é muito importante para a recuperação da hemicelulose (BRIENZO, 
2010) e tratamentos com ácidos agem severamente removendo a hemicelulose e partes 
da celulose e lignina. Outro pré-tratamento muito usado é com reagentes oxidantes, 
por exemplo os perácidos que apresentam um elevado potencial oxidante e que age na 
remoção de lignina (BRASILEIRO et. al., 2001). 

Para os tratamentos biológicos é utilizado na sua maioria enzimas lignoliticas, que 
são produzidas pelos fungos da podridão branca geralmente como metabolito secundário 
(FARIA,2011). Entre essas enzimas, que tem a capacidade de oxidar a lignina, tem-se 
as lacases que apresentam um potencial redox suficiente para a remoção de elétrons da 
estrutura fenólica, e assim aumentando sua taxa de deslignificação (BAO et. al., 1994). 
Para a extração da fração da hemicelulose e celulose é utilizada as xilanases e celulases, 
respectivamente. Ambas as enzimas fazem parte de um extrato enzimático que agem em 
sinergismo, fazendo a hidrolise da xilana em xilose, e de frações da celulose em glicose 
(TALEBNIA et al., 2010). 

Este trabalho utilizou duas metodologias com o intuito de obter-se uma maior 
deslignificação e fracionamento da biomassa lignocelulósica. A biomassa do eucalipto 
VCP (Votorantin Celulose e Papel), foi submetido ao pré-tratamento com ácido peracético, 
diferenciando-os em tempo e temperatura de tratamento. Com o eucalipto UFT, foi aplicado 
o pré-tratamento com NaOH a 5% (m/v) e posteriormente foi também tratado com lacase. A 


hidrólise das biomassas foi realizada com coquetel enzimático Cellic Ctec2º. 


21 MATERIAIS E MÉTODOS 


2.1 Materiais 


Os cavacos de eucalipto foram obtidos na indústria Votorantim celulose e papel 
(VCP - atualmente Suzano) da cidade de Jacareí/SP. E o eucalipto (Eucalyptus grandis) 
pré-tratado a 150º C foi cedido pelo laboratório de Tecnologia e Utilização de Produtos 
Florestais (LabTecFlor) do Curso de Engenharia Florestal da Universidade Federal do 
Tocantins. Os materiais foram moídos a 20 mesh e guardados em sacos plásticos em 
temperatura ambiente para a execução do experimento. 
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2.2 Tratamento de Eucalyptus VCP com ácido peracético 


O ácido peracético foi produzido a partir da mistura de acetato de etila com peróxido 
de hidrogênio na razão de 1:1 (v/v) na presença de 1% do catalisador ácido sulfúrico. 
Para o pré-tratamento, 3 g de massa seca da serragem de eucalipto foi acondicionado 
em Erlenmeyer de 250 ml com adição de 30 ml da solução de ácido peracético (PAA) 
alcançando uma concentração final de 10% (m/v). O frasco foi fechado com rolha de modo 
que fosse gerada uma pressão no seu interior. 

O pré-tratamento foi conduzido em temperaturas de 60º C e 90 ºC. Das amostras 
pré-tratadas a 60ºC foram retiradas amostras em tempos de 24, 48 e 72h para analisar o 
tempo de pré-tratamento. No pré-tratamento a 90 ºC, a reação foi encerrada com 5 horas 
na tentativa de otimizar o tempo com o aumento da temperatura. 

As amostras pré-tratadas foram filtradas em filtro de cerâmica porosa sinterizada nº 
2 à vácuo para separar a fração líquida da biomassa sólida. A biomassa sólida foi lavada 
exaustivamente com água destilada estéril até pH neutro e retirada do excesso de ácido 
peracético. As amostras foram então mantidas em temperatura de 60 “C para secagem até 


peso seco. 


2.3 Tratamento com hidróxido de sódio 


Este tratamento foi realizado com eucalipto doado pelo LabTecFlor que já havia 
sido tratado termicamente a 150 “C por uma hora. O material moído a 20 mesh (1269) foi 
acondicionado em Erlenmeyer de 1000 mL com 450 mL de NaOH a 5% (m/v) e autoclavado 
a 120 ºC, 1 atm por 1h. O material foi então neutralizado com ácido clorídrico (1,25 M) e 
filtrado sob papel de filtro Whatman nº.1 sob vácuo. O resíduo sólido foi lavado com água 
destilada e levada para estufa a 60 “C para secagem e devidamente identificada 


2.4 Tratamento com Lacase 


O extrato de Lacase, usado neste trabalho foi cedido pelo Laboratório de 
Microbiologia e Bioquímica | da Escola de Engenharia de Lorena da Universidade de São 
Paulo (EEL-USP). Em tubos de ensaio de 50 mL foi acondicionado 0,2g de massa seca 
de eucalipto pré-tratado com 5% de NaOH com 5 mL da solução tampão acetato de sódio 
(pH 5; 50 mM) e mantido a 30 “C por três dias. Então, o volume foi completado para 10 mL 
juntamente com 10 U/mL de lacase e a reação enzimática foi acompanhada por dois dias 
em banho maria a 30º C. Para comparação dos resultados, foram montados experimentos 
na presença de ABTS 30mM e o branco sem a presença de lacase e ABTS. 

Para analisar a efetividade dos tratamentos, todas as amostras de eucalipto pré- 
tratadas foram submetidas a hidrólise com enzima comercial Cellic Ctec2 e caracterizadas 
quimicamente antes e depois de cada tratamento. A equação 1 foi usada para calcular a 
recuperação da biomassa sólida. A equação 2 foi utilizada para calcular a quantidade em 
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gramas dos componentes (lignina, glicana e xilana) por 100g de eucalipto. 
R% = x 100 (1) 
PSI 


R%= Rendimento; 
PSi= Peso seco inicial; 
PSf= Peso seco final 
COpt RM 


Ep = SEA. (: 


1ÕO 


Cg= Composição em g por 100 g de eucalipto; 
R%= Rendimento; 
CQpt= Componente químico do pré-tratado; 


2.5 Caracterização química da serragem de eucalipto 


2.5.1  Quantificação de extrativos 


Em saquinhos de papel filtro Whatman nº.1 foi adicionado 1,5 g da serragem de 
eucalipto (VCP) in natura (base seca) após secagem a 60 “ºC para determinar a umidade. 
Os extrativos das amostras foram determinados depois que as amostras passaram por 
processo de extração em aparelho Soxhlet com etanol 95% por 6 horas. O teor de extrativos 
foi então determinado pela pesagem do material após secagem da serragem de eucalipto 
a 60 “C até peso constante. 


2.5.2 Análise da composição química 


As amostras de eucalipto in natura e pré-tratada, foram caracterizadas quimicamente 
pelo método descrito por Ferraz (2000). A composição química dos materiais foi realizada 
por análise em triplicata a partir de 330 mg (massa seca) de amostra em tubo de ensaio. 
Em cada tubo foram adicionados 3 mL de H,SO,a 72% (m/v) seguido de incubação em 
banho maria a 30 ºC por 1 hora. O material foi transferido para frasco Erlenmeyer de 250 
mL adicionando 79 mL de água destilada para se obter uma concentração final de 4% (m/v) 
de H,SO, e autoclavado a 121ºC, 1 atm por 1 hora. 

O material hidrolisado foi resfriado em temperatura ambiente e filtrado em filtro de 
cerâmica porosa nº 3, para a determinação de lignina insolúvel por gravimetria, e a lignina 
solúvel por leitura em espectrofotômetro a 205 nm (O = 105 L.g!.cm”) diluindo 10 vezes. 
Os volumes das soluções filtradas foram completados para 200 ml com água destilada 
em balão volumétrico. Uma alíquota de cada solução foi filtrada em filtro millipore 40 um 
e utilizada para análise das concentrações de açúcares monoméricos gerados glicose, 
xilose, arabinose e ácido acético. As análises desses compostos foram realizadas em um 
equipamento HPLC (CLAE — Cromatografia Líquida de Alta Eficiência - WATERS) equipado 
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com um detector de índice de refração e uma coluna HPX87H (BIORAD) a 35 ºC, com 
fase móvel composta por ácido sulfúrico 5 mM e um fluxo de 0,6 mL.min-1. Os açúcares 
foram expressos como seus respectivos polissacarídeos de acordo com seus fatores de 
hidrólise: 0,90 para glicana, 0,88 para xilana e arabinosil e 0,72 para acetil. As áreas dos 
cromatogramas foram convertidas em concentração (g/L) utilizando curvas de calibração 


previamente construídas com soluções padrão de glicose, xilose, arabinose e ácido acético. 


2.5.3 Determinação de cinzas 


Cadinhos de porcelana foram mantidos em muflas por 3 horas a 575 ºC, depois 
colocados em dessecador para estabilizar a temperatura antes da pesagem. Cerca de 
0,4g (massa seca) de eucalipto moído foi adicionado em cada cadinho e levado a mufl 
tampado. Uma rampa de aquecimento determinada por Sluiter (2008) foi utilizada até atingir 
575 “C em 1 hora. Os cadinhos foram destampados e aguardou-se a mufla esfriar para 
transferir os cadinhos para um dessecador para estabilizar a temperatura e em seguida 
feita a pesagem das cinzas (SLUITER, 2008). 


2.6 Hidrólise enzimática dos pré-tratamentos 


A hidrólise enzimática foi realizada com extrato comercial de enzimas Cellic Ctec2 
(Novozymes) 10 FPU/g (unidades de papel de filtro/grama de biomassa), empregando 3g 
das diversas polpas pré-tratadas. 

As amostras com 3 dias de pré-tratamento com ácido peracético a 60ºC e com 
lacase foram adicionadas cargas fixas de enzimas com 10 FPU/g diluídas em tampão 
acetato de sódio (50 mM; pH 5). As amostras pré-tratadas com ácido peracético a 90ºC 
foram avaliadas as melhores cargas enzimáticas entre 7,5, 15 e 22,5 FPU/g diluídas no 
mesmo tampão acetato de sódio. 

Os experimentos foram conduzidos em tubos de ensaios de 50 ml, contendo 10 
ml de meio reacional, sob agitação de 120 rpm e temperatura de 50ºC. Amostras foram 
retiradas periodicamente em 6, 24, 48 e 72 h e aquecias em água fervente (100ºC) por 
5 minutos para interromper a reação. As amostras foram centrifugadas a 5000 x g por 5 
minutos e os sobrenadantes foram congelados para análises conforme item 2.6.2. 

A equação 3 calcula a efetividade na extração da lignina e xilana das amostras de 
eucalipto. 

3 ] LOpt e Rb : T 
E=(1 Com | + 100 (3 


“arc! 


E= Extração 
CQpt= Componente químico do pré-tratado 
CQin= Componente químico in natura 


R%= Rendimento 
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A hidrólise enzimática da celulose e xilana foi estimada como a relação entre a 
concentração de glicana obtida após a hidrólise enzimática e a concentração de glicana 
presente no substrato inicial, de acordo com a equação 4: 


CE (a) 100 4) 
Emi | 1 (4) 


CE: Conversão enzimática 

Cg: Concentração de glicana ou xilana obtida na hidrólise 

Fc: Fator de conversão de celulose em glicose (0,9) e xilana em xilose (0,88) 
Vr: Volume da reação enzimática 

C: Porcentagem de glicana ou xilana nas amostras in natura 

M: Massa seca das amostras 


31 RESULTADOS E DISCUSSÃO 


3.1 Composição química da serragem de eucalipto in natura 


A caracterização química de ambas biomassas com relação a porcentagem de 
lignina, glicana, xilana, acetil, extrativos e cinzas são apresentadas (Tabela 1). As biomassas 
foram obtidas da indústria Votorantim Celulose e Papel (VCP) e da Universidade Federal 
do Tocantins (UFT), esta pré-tratada com 5% de NaOH, a qual não foi possível realizar a 


quantificação de extrativos em consequência de já ter sido um m terial pré-tratado. 


Componente (%) EU-VCP EU-UFT 
Lignina 30,9 + 1,0 37,9+0,4 
Glicana 41,3+0,7 47,3+0,6 
Xilana 13,2+0,4 9,1+0,2 

Arabinosil 0,3 0,2 
Acetil 3,0+0,1 0,1 
Cinzas 0,4 + 0,2 0,4 

Extrativos 2,1+0,1 0 
Total 91,2+1,2 95,4 + 0,4 


VCP — Votorantim Celulose e Papel; UFT — Universidade Federal do Tocantins 


Tabela 1. Composição química do eucalipto in natura doado pela VCP e do eucalipto (UFT) pré-tratado 
com 5% (m/v) de NaOH. 


A porcentagem de lignina da biomassa UFT foi maior que da VCP, apresentando 
37,9 e 30,9%, respectivamente. Estes valores são considerados como o total entre lignina 
solúvel e insolúvel. O inverso foi obtido no conteúdo de xilana e arabinosil, para biomassa 


VCP de 13,2 e 0,3%, enquanto a biomassa UFT foi menor em quase 5% em xilana e 
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pouca diferença no teor de arabinosil de 0,1%. Ainda com relação a xilana, que é um 
heteropolímero, o teor de xilana da biomassa UFT apresentou menor acetilação (0,1%), 
enquanto a biomassa da VCP maior acetilação foi encontrada (3,0%). A variação na 
porcentagem dos componentes das biomassas pode estar ligada nas diferentes espécies 
e linhagens de eucalipto avaliadas (Tabela 1). As porcentagens de lignina, glicana, xilana, 
grupos acetil e arabinosil se apresentam com clareza bastante diferentes. A quantidade 
de glicose nas biomassas de diferentes origens, considerada como um polímero infinit 

e denominada glicana não ficou diferente. A biomassa VCP apresentou 41,3%, enquanto 
47,3% foi encontrada na biomassa UFT, diferença que pode estar relacionada ao pré- 
tratamento com NaOH realizado antes da caracterização. 

Os dados encontrados estão em consonância com a caracterização química realizada 
por Kundu (2021), em que o eucalipto apresentou 40,51% de glicana, 22,12% de xilana, 
1,06% de arabinosil e 33,53% de lignina. A baixa acetilação da biomassa UFT pode ter sido 
ocasionada pela remoção dessa substância durante o processo de pré-tratamento como 
observado por Shimizu (2018). O autor observou que antes do pré-tratamento alcalino, a 
biomassa in natura tinha 3,42% de acetil e após o tratamento alcalino a 5% (m/v) a fração 
de acetil diminuiu para 0,86%. 


3.2 Composição química do pré-tratado com ácido peracético 


A temperatura e o tempo foram os fatores que mais influenciaram na deslignificaçã 
da biomassa VCP após o pré-tratamento com ácido peracético (PAA). Além disso, notou- 
se uma diminuição constante da recuperação de sólido com o aumento do tempo de pré- 
tratamento. Dessa maneira, é possível afirmar que o tempo influencia diretamente no 
rendimento da remoção de lignina. 

A recuperação de sólidos foi calculada pela massa remanescente após filtraçã 
da biomassa pré-tratada com PAA. A queda no rendimento está atribuída à remoção de 
xilana e lignina da biomassa analisada, sendo que a porcentagem de glicana se manteve 
praticamente constante após os pré-tratametos (Tabela 2). O cálculo da quantidade dos 
componentes da biomassa foi realizado considerando o rendimento na recuperação de 
sólidos, o total de lignina, glicana e xilana presentes no sólido recuperado e suas quantidades 
in natura antes do pré-tratamento. A lignina apresentou uma queda de 10,3% dentro de 
72h a 60ºC e a xilana de 7,8% dentro deste mesmo pré-tratamento. Em temperatura mais 
elevada (90 ºC), a deslignificação aumentou para 13,4% e a hidrólise da xilana para 9,4% 
com 5h de pré-tratamento. Isso implica que o aumento da temperatura diminui o tempo de 


processamento de degradação de lignina e xilana. 


Engenharia química: Desenvolvimento de novos processos e produtos 2 Capítulo 5 60 


Composição (%) Remoção % 


Amostras er Lignina  Glicana Xilana Lignina Glicana Xilana 
In natura 100 30,87 41,34 13,22 - - - 
24h 60 ºC PT-PAA 84,8 23,72 41,44 9,01 23,16 - 31,85 
48h 60 “C PT-PAA 78,0 22,18 42,16 5,89 28,16 - 55,4 
72h 60 ºC PT-PAA 75,6 20,61 40,88 5,41 33,25 1,12 59,03 
5h 90 ºC PT-PAA 72,6 17,50 42,62 3,78 43,32 - 71,43 


PT-PAA — Pré-tratamento com ácido peracético 


Tabela 2. Composição de lignina, xilana, glicana (%) e taxa de remoção (%) da biomassa VCP após 
pré-tratamento com PAA em diferentes condições de temperatura e tempo. 


Não foi observado aumento expressivo no conteúdo de glicana em todos os pré- 
tratamentos, mas na remoção de lignina e xilana foi bastante expressiva. A remoção de 
lignina e xilana alcançou a porcentagem de 43,32% e 71,43% a 90 ºC em 5h de tratamento 
com PAA, enquanto na temperatura de 60 ºC por um período de 72h alcançou 33,25% 
para lignina e de 59,03% para xilana, observando mais uma vez que em temperatura mais 
elevada se alcança maiores porcentagens de remoção destes componentes em menor 
tempo de tratamento. 

Aumento na porcentagem de glicana em madeira dura pré-tratada com ácido 
peracético foi observado por Kundu et al. (2021), o que foi atribuído a perda de lignina e 
xilana no pré-tratamento. Neste pré-tratamento observou-se uma queda do rendimento de 
massa das amostras de eucalipto, com o decorrer do tempo e aumento da temperatura, 
que segundo Kundu et. al. (2021) essa perda de massa pode ser atribuída a remoção de 
lignina com contribuições do açúcar hemicelulósico. Os fatores observados nesse estudo, 
não foram diferentes dos fatores encontrados no experimento realizado por Kundu et. al., 
(2021), do qual tanto para madeira dura (eucalipto) quanto para madeira macia (pinus), 
houve uma indicação clara de remoção de lignina em função do tempo do pré-tratamento. 
No experimento realizado por Kundu et. al. (2021) observou-se que o teor de glicana 
permaneceu o mesmo em até 5h de pré-tratamento, mas cerca de 75% da xilana e 90% da 


lignina foram removidas da biomassa sólida. 


3.3 Extração de xilose e glicose de VCP pré-tratada com PAA a 60 “C 


A hidrólise enzimática de glicana e xilana das amostras de eucalipto VCP pré- 
tratadas com ácido peracético em diferentes tempos a 60 “C foi conduzida com o coquetel 
enzimático comercial Cellic Ctec2 com atividade de 10 FPU/g de eucalipto recuperado por 
até 72h, 120 rpm a 50 ºC (Figuras 1 e 2). Os dados foram obtidos por análise cromatográfic 
CLAE (Cromatografia Líquida de Alta Eficiência) conforme metodologia descrita em 2.6.2 

A remoção de xilana da biomassa de VCP apresentou um aumento com o tempo 
de pré-tratamento a 60 “C com ácido peracético (Tabela 2). Independente desse aumento, 
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a hidrólise enzimática não apresentou uma grande variação na quantidade de xilose 
extraída após 24h de tratamento, sendo que pequeno aumento de monômeros de xilose 
foi observado no pré-tratamento com ácido peracético por 1 dia dentro de 48h de hidrólise, 
alcançando 20% de extração (Figura 1). 


Extração xilose PT-PAA 6050 
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Figura 1. Hidrólise enzimática de xilana após pré-tratamento em diferentes tempos com ácido 
peracético a 60 ºC. Hidrólise com coquetel enzimático Cellic Ctec2 a 50 ºC, 120 rpm. 


Processo semelhante ocorreu com a extração de glicose, em que foram obtidos 7% 
deste monômero da biomassa VCP recuperada após pré-tratamento com ácido peracético 
por 1 dia (Figura 2). O processo de hidrólise com as biomassas pré-tatadas por 2 e 3 


dias não apresentou diferença na extração deste monômero em que 4% de glicose foram 
obtidas. 


Extração de glicose PT-RAA 6090 
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Figura 2. Hidrólise enzimática de xilana após pré-tratamento em diferentes tempos com ácido 
peracético a 60 ºC. Hidrólise com coquetel enzimático Cellic Ctec2 a 50 ºC, 120 rpm. 
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A remoção de constituintes da parede celular como lignina e xilana ocasiona um 
aumento da porosidade da parede celular aumentando a acessibilidade das enzimas 
e, com isso, um maior rendimento da sacarificação (SINGH, 2017). A sacarificação das 
amostras após PT-PAA a 60ºC não foi eficiente tanto quanto o que foi observado por Wen 
et. al. (2021) realizando o mesmo tipo de extração com biomassa de álamo que dentro de 


2h recuperaram 90,1% de glicose e 84,9% de xilose nesta hidrólise enzimática. 


3.4 Extração de xilose e glicose de VCP pré-tratada com PAA a 90 “C 


O aumento da temperatura no pré-tratamento da biomassa VCP parece ter 
contribuído com o processo de hidrólise de xilana e glicana (Figuras 3 e 4). Neste processo 
de hidrólise foi utilizada várias cargas da enzima comercial Cellic Ctec2 na razão de 7,5, 
15 e 22,5 FPU/g para verificar a eficiência de extração de xilose e glicose, uma vez que 
a remoção de lignina e hemicelulose foi maior em menor tempo de tratamento (Tabela 2). 

A cinética enzimática de hidrólise de xilana aumentou com a carga de Cellic Ctec2 
alcançando 40% de extração de xilose dentro de 48h usando 22,5 FPU/g de biomassa 
remanescente após tratamento com PAA a 90 “C por 5h (Figura 3). No mesmo tempo de 
hidrólise, as cargas de enzimas de 7,5 e 15 FPU/g extraíram xilose abaixo de 30%. Com o 
aumento do tempo de hidrólise para 72h, nenhum efeito maior de extração foi observado. 
Os resultados dentro desta condição de pré-tratamento se apresentaram muito mais 
promissor para o processo de recuperação de xilose e glicose. 
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Figura 3. Hidrólise enzimática de xilana com Cellic Ctec2 com diferente FPU (unidades de papel de 
filtro) após pré-tratamento em diferentes tempos com ácido peracético a 90 “C por 5h. Hidrólise com 
coquetel enzimático Cellic Ctec2 a 50 ºC, 120 rpm. 
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Não obstante, a extração de glicose após esse tratamento com PAA a 90 “C por 5h 
foi bem mais expressiva, chegando a 15% em 48h com coquetel enzimático de 22,5 FPU/g 
de VCP (Figura 4). O dobro de extração foi observado quando comparado com a biomassa 
VCP pré-tratada a 60 “C por até 3 dias (Figura 2), sendo que a metade do valor de extração 
de glicose foi obtida. 
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Figura 4. Hidrólise enzimática de glicana com Cellic Ctec2 com diferente FPU (unidades de papel de 
filtro) após pré-tratamento em diferentes tempos com ácido peracético a 90 ºC por 5h.Hidrólise com 
coquetel enzimático Cellic Ctec2 a 50 ºC, 120 rpm. 


Cellic Ctec2 º é uma mistura de enzimas hemicelulolíticas composta por xilanase, 
endoglicanase, exoglicanases e B-glicosidase que atuam na hidrólise de xilana e celulose 
(glicana). A dosagem de celulase de 10-30 (FPU/g de biomassa) é bastante usado em 
pesquisas de extração de açúcares monoméricos, resultando em hidrólise eficiente com 
rendimento de glicose alto, levando em consideração um tempo razoável de hidrólise de até 
72heo custo da enzima (TALEBNIA etal., 2010). Aatuação enzimática é mais efetiva quando 
a extração da lignina é mais eficiente nos processos de pré-tratamento. O pré-tratamento 
é primordial para inicialmente remover lignina e hemicelulose da matriz lignocelulósica, 
diminuindo a cristalinidade da celulose e aumentando a área superficial e porosidade da 
biomassa. Consequentemente, pode também liberar alguns açúcares fermentescíveis por 
despolimerização da fração hemicelulósica (ZABED et al., 2017). Sendo assim, o melhor 
pré-tratamento é aquele que deslignifica o máximo possível procurando manter a fração 
hemicelulósica praticamente intacta. 

Serragem de E. grandis foi submetida a hidrólise por uma preparação comercial 
de celulases de Trichoderma reesei após pré-tratamento biológico com vários fungos por 


Castoldi et al., 2014. A hidrólise foi grandemente aumentada na presença de celulase de 
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T.reesei e que na maioria dos casos a cinética enzimática não cessou em 48h e até mesmo 
acelerou este tempo. Preparações comerciais de celulases e xilanases diferenciam-se na 
capacidade de hidrólise de biomassa pré-tratada com soluções ácidas e/ou alcalinas. De 
fato, Accellerase*1500 foi mais efetivo, duas vezes mais, que AccelleraseXC na extração 
de glicose e xilose de fibras de frutas de óleo de palma pré-tratada com ácido peracético 
alcalino a 20 “€ com NaOH 4% m/v por 12h (PALAMAE et al., 2017). 


3.5 Pré-tratamento com NaQOH, lacase e hidrólise com Cellic Ctec2 


O eucalipto UFT pré-tratado com hidróxido de sódio a 5% (m/v) e autoclavado a 120 
ºC, 1 atm por 1h foi submetido a outro pré-tratamento com lacase (10 U.mL'.9'!) a 30 ºC por 
2 dias antes da adição do coquetel enzimático Cellic Ctec2. Para comparação, o mesmo 
experimento foi realizado na presença de ABTS (30 mM) e sem a presença de lacase e 
ABTS foi realizada experimento de hidrólise somente com Cellic Ctec2 como referência 
(Figuras 5 e 6). Os extratos obtidos da filtração da biomassa UFT tratada com lacase e 
ABTS foram analisados por método de varredura em diferentes comprimentos de onda no 
UV/Vis e apresentaram um pico em 215 nm em igual intensidade de absorbância. Lacases 
são enzimas oxidativas e degradam lignina e compostos fenólicos sem ou com a presença 
do mediador ABTS. De qualquer forma, os espectros após os tratamentos com lacase e 
lacase com ABTS foram muito semelhantes com relação a absorbância. 

Dentro destes tratamentos, as hidrólises de xilana apresentou uma extração de 
xilose mais efetiva quando a biomassa UFT foi submetida a deslignificação com lacase 
na presença de ABTS. Após esse tratamento, o uso do coquetel enzimático de celulases 
favoreceu maior extração de xilose, alcançando 37,1% em comparação aos experimentos 
somente com lacase (30,6%) e somente com Cellic Ctec2 (20,9%). Cineticamente, o tempo 
de 72h de hidrólise foi o de maior rendimento (Figura 5). 

A hidrólise de celulose em monômeros de glicose foi muito semelhante em todos os 
tratamentos (Figura 6). O tratamento com lacase na presença de ABTS apresentou uma 
extração de 10%, enquanto 9% foi alcançado somente com lacase. A extração de glicose 
com Cellic Ctec2 sem o pré-tratamento com lacase e/ou ABTS chegou em 6%. Analisando 
os resultados anteriores, da biomassa VCP pré-tratada com ácido peracético a 90 “C, a 
extração de glicose superou em mais de duas vezes os tratamentos da biomassa UFT com 
NaOH 5% (figura 4) 


Engenharia química: Desenvolvimento de novos processos e produtos 2 Capítulo 5 RE 


Extração de xilose 


ao 
ni 
30 E ai 
ch ai 
= 0 a i==— & E 
E e pa e. 
E Ep 
En “E- 
"á 
od 
o 1z | E) ER &0 TE 
Tempo [b) 
—4— LaHABTS+ tec Cimo +— Lac 


Figura 5. Extração de xilose da biomassa UFT pré-tratada com NaOH 5% (m/v) 


Extração de glicose 


15 
É io ——A 
Es a me —M 
ES E e no 
ú 
ú 12 2 36 <8 so 7. 
Tempo (b) 
—— LacaBTs-Cieci Creci +— Lac 


Figura 6. Extração de glicose da biomassa UFT pré-tratada com NaOH 5% (m/v) 


Muitos pré-tratamentos de matérias-primas lignocelulósicas são baseadas 
em processos alcalinos. A utilização de soluções diluídas de hidróxido de potássio, 
cálcio e amônio são agentes alcalinos disponíveis para o pré-tratamento de biomassa 
lignocelulósica, mas o hidróxido de sódio é o mais estudado. O pré-tratamento alcalino 
pode grandemente melhorar a digestibilidade de celulose e a degradação de açúcares é 
menor que no tratamento ácido (TALEBNIA et al., 2010). 
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Os pré-tratamentos apresentaram diferenças em procedimentos com soluções 
alcalinas e ácidas. O ácido peracético deslignificou a biomassa lignocelulósica VCP 


grandemente com o aumento da temperatura com consequente diminuição de tempo de 
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exposição ao ácido. A remoção de lignina foi maior dentro destas condições, mas diminuiu 
o rendimento da recuperação da biomassa após o pré-tratamento ácido. 

A hidrólise da biomassa com Cellic Ctec2, após o pré-tratamento ácido em 
temperatura mais elevada e menor tempo, foi a que promoveu maior extração de glicose e 
xilose. Dentro destes monômeros de açúcares, a xilose foi hidrolisada em parte com o pré- 
tratamento ácido, mas a glicose permaneceu a mesma em composição antes da hidrólise 
enzimática. 

A hidrólise enzimática da biomassa UFT; eucalipto pré-tratado por autoclavagem 
a 121 ºC por 1h na presença de hidróxido de sódio; não apresentou o mesmo resultado 
do pré-tratamento ácido em pouco tempo de cinética, somente alcançando os resultados 
cinéticos enzimáticos em longo tempo. 
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RESUMO: A preocupação com a eficiência 
energética e o meio ambiente é cada vez 
maior, pois sabe-se que quanto mais energia se 
consome, maior a emissão de gases de efeito 
estufa, que aumentam a temperatura média de 
nosso planeta. O objetivo deste artigo é avaliar 
a viabilidade tecnicoeconômica do uso do calor 
residual de chillers alimentando um chiller tipo 
bomba de calor para a geração de água quente 
para uso em hospitais, que são consumidores 
intensivos de energia elétrica, utilizada na geração 
de água gelada, e combustível fóssil na geração 
de água quente, e, com isso, eliminar ou reduzir 
substancialmente o consumo de gás liquefeito de 
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petróleo (GLP) ou gás natural nos aquecedores 
de água. Para tanto, foi avaliado o perfil de carga 
térmica de um sistema de ar condicionado e o 
perfil de uso de água quente de um hospital- 
escola em Campinas, foi realizada a comparação 
da simultaneidade dos perfis de carga e o uso de 
software de modelagem energética, que mostra 
a redução do acionamento dos aquecedores de 
água que queimam gás liquefeito de petróleo. 
A partir dos resultados dessa modelagem 
energética foi obtido o novo consumo de 
combustível fóssil do hospital nos aquecedores 
de água e o novo consumo de energia elétrica 
no conjunto de chillers existentes e com o chiller 
bomba de calor. A partir daí foi feita a análise 
de viabilidade econômica da implantação desse 
sistema de recuperação de calor residual dos 
chillers existentes e a redução da emissão anual 
de COZ. 

PALAVRAS-CHAVE: Calor Residual; Chillers; 
Bomba de Calor; Economia de Energia; 
Recuperação de Calor. 


TECHNICAL AND ECONOMIC 
FEASIBILITY OF USING RESIDUAL 
HEAT FROM CHILLERS FOR HEATING 
PROCESS IN HOSPITALS 


ABSTRACT: The concern with energy efficiency 
and the environment is increasing, as it is known 
that the more energy is consumed, the greater 
the greater is the emission of greenhouse gases, 
which increase the average temperature our 
planet. The target of this article is to evaluate the 
technical-economic feasibility of using waste heat 
from chillers feeding a heat pump type chiller for 
the generation of hot water for use in hospitals, 
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which are intensive consumers of electricity for the generation of cold water and fuel fossil to 
generate hot water and thereby eliminate or substantially reduce the consumption of liquefied 
petroleum gas (LPG) or natural gas in water heaters. Therefore, the thermal load profile of 
an air conditioning system and the hot water use profile of a teaching hospital in Campinas 
were evaluated, the comparison of the simultaneity of the load profiles and the use of energy 
modeling software that show the reduction in the activation of water heaters that burn liquefied 
petroleum gas. From the results of this energy modeling, the new fossil fuel consumption 
of the hospital in the water heaters and the new electrical energy consumption in the set of 
existing chillers and with the heat pump was obtained. From there, the economic feasibility 
analysis of the implementation of this waste heat recovery system from the existing chillers 
and, the reduction of the annual CO2 emission was carried out. 

KEYWORDS: Waste Heat; Chillers; Heat Pump; Energy Saving; Heat Recovery. 


11 INTRODUÇÃO 


No mundo inteiro, a preocupação com a eficiência energética e o meio ambiente é 
cada vez maior, pois sabe-se que quanto mais combustíveis fósseis se consome, maior a 
emissão de gases de efeito estufa. Os sistemas de ar condicionado e de aquecimento de 
água são dois grandes vetores de consumo de energia em edifícios comerciais — são os 
principais vetores energéticos do edifício, principalmente em edifícios corporativos para 
uso hospitalar. 

Nesse ambiente, a operação funciona 24 horas, durante sete dias da semana, e 
os sistemas de ar condicionado e de aquecimento de água precisam ser mantidos ligados 
durante esse período. Em sua maioria, os hospitais utilizam sistemas de ar condicionado 
através de uma central de água gelada (CAG) para a climatização de ambientes classificados 
e sistemas de aquecedores com queima de combustível fóssil para aquecimento de água 
para banho e cozinha, principalmente. Observa-se, portanto, o uso simultâneo dos sistemas 
de resfriamento e aquecimento. 

Aí surge um questionamento: o rejeito térmico do condensador de um chiller, que é 
expurgado para o meio ambiente ou através de uma torre de resfriamento, não pode ser 
reaproveitado como calor potencial para uso em outro sistema? 

Estudos mostram o desenvolvimento de sistemas de geração de água gelada e 
água quente integrados, com o aproveitamento de parte ou todo o rejeito térmico das 
unidades resfriadoras de líquido de uma central de água gelada, e a eliminação total ou 
parcial do uso de combustível fóssil, como gás natural ou óleo combustível, de caldeiras ou 
aquecedores de passagem, melhorando-se a eficiência energética de um hospital. Assim, 
esta pesquisa tem como objetivo avaliar a viabilidade técnica e econômica da implantação 
de um sistema de um chiller bomba de calor para aquecimento de água para uso em 
um hospital. Esse procedimento é feito através do aproveitamento do rejeito térmico do 
chiller da central de água gelada do hospital, ao invés de gerar essa água quente através 
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de sistemas de aquecimento que queimam combustível fóssil, ou seja, gás liquefeito de 
petróleo (GLP), reduzindo-se a emissão de CO2 para a atmosfera e, consequentemente, o 
global warm potential (GWP). Essa avaliação é feita a partir do levantamento dos perfis de 
carga térmica de ar condicionado e de aquecimento desse hospital. 


21 REVISÃO BIBLIOGRÁFICA 


Esta pesquisa buscou entender essa situação em edifícios hospitalares, tanto em 
termos de consumo de energia elétrica, quanto do consumo de combustíveis. A pesquisa 
foi estendida até a China, país com a maior população do mundo e o que mais cresce 
anualmente em termos econômicos. De acordo com Chiang et al. (2017), recentemente 
a China se tornou o maior consumidor de energia do mundo, e o uso de energia em um 
hospital é duas vezes maior que o consumo de um edifício público. 

Esse alto consumo de energia nos edifícios hospitalares ocorre devido às cargas 
de aquecimento e resfriamento, que acontecem de forma contínua, principalmente com a 
climatização dos ambientes e o aquecimento de água sanitária. Os edifícios dos hospitais 
usam uma variedade de fontes de energia para aquecimento e resfriamento, incluindo 
eletricidade, energia renovável, gás natural e combustível fóssil. Dentre eles, o sistema de 
climatização é o que mais utiliza energia elétrica, respondendo por 50% de todo o consumo, 
conforme a Figura 1. 
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Figura 1 — Consumo de energia nas utilidades de um hospital típico 


Fonte: Chiang et al., 2017, p. 2 


Já os combustíveis fósseis são mais comumente usados como fonte de 
abastecimento no sistema doméstico de aquecimento de água. De acordo com Chiang et 
al. (2017), mostrado na Figura 2, o uso de lavar e secar contribui com 58% do consumo 
total de combustível; aquecimento e uso de cozinha contribuem com 23%; e o uso de 


chuveiro contribui com 17% do consumo total de combustível. 
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Figura 2 — Consumo de combustível fóssil nas utilidades de um hospital típico 

Fonte: Chiang et al., 2017, p. 2 


E o que fazer com o rejeito térmico do condensador de um chiller da central de 
geração de água gelada do sistema de ar condicionado? Esse calor é simplesmente 
expurgado para o meio ambiente através de uma torre de resfriamento ou diretamente para 
o ar? 

De acordo com Jia (2006), existem dois métodos de recuperação de calor de um 
chiller. um método utiliza um trocador adicional, fora do condensador resfriado a água, 
e em vez de rejeitar o calor para a torre de resfriamento, este é recuperado da água do 
condensador, conforme Figura 3; no outro método, a recuperação de calor pode ocorrer 
diretamente do fluido refrigerante quente no circuito de refrigeração usando um feixe de 
condensador adicional, conectado em paralelo com o condensador padrão, conforme 
Figura 4. Deve-se lembrar que um dos requisitos para a recuperação de calor em um 
sistema de produção de água gelada é haver aquecimento e o resfriamento simultâneos. 


Figura 3 — Recuperação de calor com trocador de calor 
Fonte: Jia, 2006 


Engenharia química: Desenvolvimento de novos processos e produtos 2 Capítulo 6 


Figura 4 — Recuperação de calor com condensador adicional 


Fonte: Jia, 2006 


Outro ponto a comentar de nossa pesquisa é que nos deparamos por diversas vezes 
com os termos bomba de calor e recuperador de calor, usados alternadamente. De acordo 
com Jia (2006), quando a função principal de um resfriador é fornecer resfriamento e uma 
parte do calor rejeitado é usada para satisfazer as cargas de aquecimento, é conhecido 
como chiller com recuperação de calor. Entretanto, quando um resfriador de líquido é 
usado para fornecer calor de seu condensador como a sua função principal, o resfriador é 
considerado uma bomba de calor. 

O resfriador foco deste artigo é a bomba de calor, e o calor disponibilizado é 
empregado no aquecimento da água, que tem utilização em aplicações comerciais e 
fornece água quente para aquecimento na faixa de 40,6ºC a 45ºC. 

Esses chillers trabalham com uma temperatura de condensação mais elevada, 
como vemos em Jia (2006) e o apresentado na Figura 5. Já saem da fábrica com recursos 
de controle de capacidade que usam a temperatura de saída da água do condensador para 
determinar qual estágio de capacidade do compressor é necessário para manter o ponto 
de ajuste da temperatura de água quente. O chiller ajustará a capacidade do compressor 


para garantir a manutenção da temperatura da água quente. 
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Figura 5 — Diferença de pressão do refrigerante 


Fonte: Jia, 2006 


Conforme Trane (2007), o controle de capacidade dos chillers com recuperação de 
calor dependerá da forma de descarregamento do resfriador. Ele pode ocorrer pelo controle 
de temperatura da água quente na saída do condensador, conforme indicado na Figura 6, 
ou na entrada do condensador, conforme indicado na Figura 7. Para compressores tipo 
scroll ou parafuso, o ideal é o controle da temperatura na saída do condensador. Essa 
estratégia procura carregar o resfriador ao máximo possível e manter a temperatura 
desejada da saída de água do resfriador, permitindo, assim, que a função principal do 
resfriador seja satisfazer a carga de aquecimento. O mesmo tipo de controle pode ser feito 
para compressores tipo centrífugo, desde que seja possível evitar a condição de surge. 


E, 


a E ma 


Figura 6 — Controle de temperatura de saída da água 


Fonte: Trane, 2007, p. 5 
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Figura 7 — Controle de temperatura de saída da água quente 


Fonte: Trane, 2007, p. 


Assim, conforme apresentado na Figura 7, o controlador de temperatura da água do 
condensador (T2) modula a válvula de desvio no laço da torre de resfriamento (V2) para 
manter uma temperatura constante de água quente. O controlador de temperatura da água 
do condensador (T2) modula a válvula de desvio no laço da torre de resfriamento (V2) para 
manter uma temperatura constante na temperatura de entrada de água no condensador. 
Este controlador também pode controlar o ventilador da torre para manter a temperatura de 
entrada da água no condensador. 

Também se faz necessário o controle de temperatura da água de saída do 
condensador (T4). O seu ajuste é feito através do controle da saída da água gelada no 
evaporador, ou seja, descarrega-se ou, carrega-se o resfriador caso seja necessário 
aumentar ou diminuir T4. Por fim a terceira malha de controle a ser inserida é para o 
controle das cargas de aquecimento através da válvula V1 controlada a partir de T1. Assim, 
uma fonte de calor auxiliar pode ser necessária se a recuperação de calor não for alta o 
suficiente para manter a temperatura de saída do condensador. 

De acordo com a Figura 15 na pesquisa de Chiang et al. (2017), um aquecedor de 
água reserva deve ser integrado ao sistema para garantir que uma fonte controlada de 
água quente esteja disponível durante os momentos em que o resfriador de recuperação 


de calor não estiver operando, como indicado na Figura 8. 
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Figura 8 — Unifilar de uma bomba de calor (WHP) integrada a uma planta de água gelada e fan & coils 
(AHU) 


Fonte: Chiang et al., 2017, p. 16 


Outra preocupação que notamos em nossa pesquisa é a questão de como configura 
o sistema de recuperação de calor. De acordo com a pesquisa Waterside Heat Recovery 
in HVAC Systems (SCHWEDLER; BRUNSVOLD, 2011), um resfriador com tubulação em 
uma posição de fluxo lateral, conforme mostrado na figura abaixo, pode ser carregado em 
qualquer condição de carga, variando seu ponto de ajuste de água resfriada. Isso ocorre 
porque o resfriador de fluxo lateral, ou sidestream, pré-resfria a água, e estes resfriadores 
têm uma temperatura de entrada de água inferior à temperatura da água das cargas de 
resfriamento. 

A principal vantagem desse arranjo é que o resfriador não precisa produzir a 
temperatura da água de abastecimento do sistema principal. Ele pode produzir qualquer 
temperatura de água necessária para atender exatamente à carga de aquecimento exigida. 
Isso permite que o resfriamento seja produzido em um nível mais alto de temperatura de 
água gelada, conforme apresentado na Figura 9. 
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Figura 9 — Arranjo tipo sidestream 


Fonte: Trane, 2007, p. 7 


Assim, de acordo com Chiang et al. (2017), um estudo feito em um dos grandes 
centros médicos públicos no Sul de Taiwan, cnamado Hospital K, mostrou que a instalação 
de um sistema de resfriador de água gelada com bomba de calor (WSHP) e integração com 
o sistema de ar condicionado existente torna-se uma solução tecnicoeconômica viável, 
conforme indicado na Figura 9. Os resultados desse retrofit mostraram uma economia 
significativa de energia e redução dos custos operacionais com um retorno financeiro do 
investimento da ordem de 1,2 anos. 


31 MATERIAIS E MÉTODOS 


Para a realização deste trabalho, estudou-se a operação atual da central de utilidades 
do Hospital de Clínicas da Universidade Estadual de Campinas (Unicamp) para a produção 
de água gelada para o sistema de ar condicionado central e produção de água quente para 
utilização em banhos, cozinha e lavanderia, e os respectivos gastos com energia elétrica 
e combustível. 

A central de produção de água gelada do hospital é composta atualmente por dois 
resfriadores de líquido tipo chiller com condensação à água e compressores tipo parafuso, 
com capacidade nominal de 300TR cada, sendo um reserva. Já a central de produção de 
água quente é composta por oito aquecedores de passagem através de queima de gás 
liquefeito de petróleo (GLP), com capacidade total de 800kW. 

Nesta avaliação, utilizou-se o perfil energético do hospital da Unicamp fornecido 
pelo hospital em termos de demandas operacionais de eletricidade, refrigeração e 
aquecimento, que estão indicados no trabalho de conclusão de curso Avaliação energética 
de implantação de um sistema de trigeração com termoacumulação mediante um estudo de 
caso (OLIVEIRA, 2017), conforme indicado no Quadro 1 abaixo. Neste trabalho, pudemos 
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observar que, ao longo de toda a operação anual, existe uma simultaneidade de consumo 


de água gelada e de água quente. 


Carga térmica Carga térmica 
Mês de resfriamento de aquecimento 
(kWh) (kWh) 

Janeiro 564.390 456.127 
Fevereiro 509.840 412.022 
Março 568.630 458.397 
Abril 361.600 378.360 
Maio 375.120 394.575 
Junho 363.800 383.765 
Julho 255.600 336.902 
Agosto 256.740 332.968 
Setembro 247.600 325.188 
Outubro 446.290 362.904 
Novembro 431.400 351.606 
Dezembro 442.890 365.718 


Quadro 1 — Perfis de carga térmica anua 


Fonte: Elaborado pelos autores 


A partir desses perfis de carga térmica e das características de capacidade e 


eficiência dos equipamentos das centrais de água gelada e água quente, foi feita uma 


simulação computacional para avaliação dos consumos de energia elétrica e térmica das 


referidas centrais térmicas. Para tanto, utilizou-se o software Trace 700, disponibilizado 


pela Trane Technologies. 
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Perfis anuais de Carga térmica 
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Figura 10 — Perfis anuais de carga térmic 


Fonte: Elaborado pelos autores 


Foi simulada a configuração das centrais térmicas atuais e, na sequência, foi 
simulada, com os mesmos perfis de carga térmica, uma central térmica de água gelada 
e água quente, em acordo com o Quadro 1 e a Figura 10, com a introdução de um chiller 
tipo bomba de calor resfriado à água, posicionado na central de água gelada atual em um 
arranjo tipo sidestream, conforme a Figura 9. Para a produção de água quente do hospital, 
considerou-se a operação do atual banco de aquecedores de passagem, que permanece 
em standby, e atuará apenas em caso de necessidade de complemento de carga, quando 
não houver produção de água gelada suficiente para a recuperação do rejeito do calor para 
a produção de água quente, conforme mostrado na Figura 9. 

Os dados de entrada solicitados pelo software Trace 700 estão indicados abaixo: 


1. Chillers da central de água gelada 


Coefficient of performance (COP): 0,683 


[º) 


. Aquecedores de passagem 


Eficiência: 83 
3. Chiller bomba de calor: 220TR de capacidade nominal 
3.1 Para a produção de água gelada 
* Capacidade de produção: 585,9kW 
* COP: 1,15 
3.2 Para a produção de água quente 


* Capacidade de produção: 596,5kW 
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* COP: 1,27 


No arranjo hidráulico do chiller bomba de calor foi utilizada a unidade com 
condensador único, conforme consta na Figura 3. Nas premissas deste estudo também 
foi considerado que a função primária desse resfriador de líquido é satisfazer a carga de 
aquecimento, ficando assim a produção de água gelada como um subproduto. Dessa forma, 
o controle de setpoint do chiller bomba de calor visa o controle de temperatura constante na 
saída do condensador, de acordo com a Figura 6. 

Os resultados dessas simulações, com seus respectivos consumos de energia 
elétrica e térmica, estão indicados no Quadro 2 a seguir, categorizados por equipamento 


das centrais de água gelada e central de água quente (CAQ). 


Consumo mensal de insumos Consumo mensal de insumos 
CAG/CAQ atuais CAG/CAQ previstos 
Mês Chillers Boilers Torres Chillers Boilers Torres Chiller Heat 
(kwh) (kWh,t) (kWh) (kWh) (kwWh,t) (kWh) Pump (kWh) 
Janeiro 70.144 679.447 9.949 7.078 2.008 6.645 135.018 
Fevereiro 6.024 613.748 8.986 6.663 1.808 6.006 121.961 
Março 70.248 682.727 9.949 7.433 1.606 6.928 135.616 
Abril 52.050 490.515 9.628 - 26.238 - 94.072 
Maio 47,25 507.101 9.869 - 27.798 - 97.547 
Junho 46.484 490.869 9.489 - 27.261 - 94.581 
Julho 30.532 505.873 8.225 - 99.546 - 68.303 
Agosto 31.509 499.964 9.350 - 125.306 - 68.588 
Setembro 30.161 488.284 8.820 - 126.141 - 66.161 
Outubro 62.336 544.917 9.949 - 5.953 - 113.522 
Novembro | 60.097 527.952 9.628 - 5.903 - 109.753 
Dezembro | 62.607 549.143 9.949 - 6.941 - 112.788 
Total 628.817 | 6.580.541 | 113.791 | 21.174 456.509 19.580 1.217.910 


Quadro 2 — Consumo de energia elétrica e térmica das centrais de utilidades 


Fonte: Elaborado pelos autores 


A próxima etapa da metodologia deste trabalho foi avaliar a viabilidade econômica 
da aplicação do resfriador tipo bomba de calor em um arranjo hidráulico tipo sidestream 
para o aquecimento da água do hospital. Foi considerado o uso complementar da carga 


de aquecimento pelos aquecedores atuais quando não houver produção de água gelada 
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suficiente para a recuperação do rejeito do calor para a produç o de água quente. 

Para tanto, buscou-se junto ao hospital as tarifas de energia elétrica e de GLP 
atualmente pagas, bem como de água, pois, conforme mencionado anteriormente, os chillers 
têm seus condensadores resfriados por água, que é resfriada em torres de resfriamento, 
o mesmo ocorrendo com o chiller tipo bomba de calor utilizado nesta simulação, conforme 
o Quadro 3. 

Os resultados dos custos operacionais com essas utilidades — ou seja, energia 
elétrica, GLP e água para reposição nas bacias das torres de resfriamento — estão 
mostradas no Quadro 4, para as duas situações, ou seja, na central de utilidades atual e 


na solução proposta. 


Tópico Valor 
Tarifa (A4-Azul)-ACR-R$/kWh 0,64 
Tarifa de GLP-(R$/kg) 3,40 
Tarifa de água-(R$/Ig) 25,00 
PCI GLP (kcal/kg) 4.924 
Carga térmica instalada de frio (TR) 300 
Carga térmica nd de aquecimento 78.700 
Tempo de operacão anual (h) 8.760 


Carga térmica anual de resfriamento (TR-h) 1.378.257 


Carga térmica anual de aquecimento 


(KWh) 4.558.031 


Quadro 3 — Premissas para a avaliação de custos operacionais (data base junho/2021) 


Fonte: Elaborado pelos autores 


No Quadro 4 estão indicados os resultados apresentados pela simulação com 
relação aos gastos operacionais com energia elétrica e energia térmica no arranjo das 
atuais centrais térmicas e com a utilização de um chiller que recupera o calor retirado 
pela água gelada utilizada nos condicionadores de ar tipo fan & coil. No Quadro 5 estão 
indicados os valores de emissão de CO2 para a atmosfera, calculados para as situações 
atual e proposta, lembrando que os valores de emissão devido ao consumo de energia 
elétrica foram calculados de acordo com o fator médio anual de t-CO2/mWh informado pelo 
Ministério da Ciência e Tecnologia, conforme valores despachados no ano de 2021. 
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Consumo anual de energia 


Insumo Situação atual Situação proposta Variação (%) 
Eletricidade (kWh) 1.526.843 2.186.186 43% 
GLP (kg) 427.859 10.005 -98% 
Água (m?) 8.741 3.645 -58% 


Gasto anual de energia 


Insumo Situação atual Custo anual (R$) Variação (%) 
Eletricidade (R$/ano) 983.287 1.407.904 43% 
GLP (R$/ano) 1.454.721 34.017 -98% 
Água (R$/ano) 218.525 91.125 -58% 
Total (R$/ano) 2.656.533 1.533.046 -42% 


Quadro 4 — Consumo de energia e custos operacionais (data base junho/2021) 


Fonte: Elaborado pelos autores 


Emissão anual de CO2 


LE Situação miual | Situação proposta [Naninção (2a) 

Eletricidade (Mwh) 193 276 43% 
GLP (ton) 1,292 30 -98% 
Total (R$/ano) 1,485 307 -79% 


Quadro 5 — Emissão de t-CO2 para a atmosfera (data base janeiro/2022) 


Fonte: Elaborado pelos autores 


Na última etapa da nossa metodologia, avaliamos o valor do capital a ser investido 
para um chiller tipo bomba de calor e demais acessórios. As despesas de capital incluíram 
a aquisição do chiller bomba de calor, do sistema de automação e a montagem do novo 
anel hidráulico sidestream, incluindo-se as bombas de água gelada e de água quente do 
condensador, lembrando que 60% desse montante foi cotado em dólar americano, com 
cotação de R$5.70/US$ (data base 07/01/2022). Na Tabela 1 e Quadro 6 estão indicados o 
fluxo de caixa e resultados financeiros desse investiment 


FLUXO DE CAIXA DO PROJETO (milhares de reais) 


Investimento Inicial (mil R$) 2350 | | 


Economia do Projeto 1123 | 1146 1,169 | 1,192 1216 1,240 1,265 1,291 1316 [139 1,370 1,397 1,425 1453 1,482 


TOTAL -2,350, 1,123 | 1,146 | 1,169 | 1,192 1216 1,240 1,265 1,291 1316 1,343 1,370 1,397 1,425 1,453 1,482 
SALDO ACUMULADO -2,350 -1,227 [a 1,088 | 2,281 3497 4,737 6,002 7293 8,609 9,952 11,321 | 12,718 | 14,143 | 15,597 | 17,079 
Py pisa do Toma 


Tabela 1 — Fluxo de caixa do projeto, vida útil de 15 anos 


Fonte: Elaborado pelos autores 
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INDICADORES FINANCEIROS DO INVESTIMENTO 


VPL fl R$) 
ROL 
THA Projeto 


Fagbiack simples + anos 


Quadro 6 — Resultados financeiros do fluxo de caixa/investimento no proje 


Fonte: Elaborado pelos autores 


A partir daí, chegou-se a um valor de capital da ordem de R$2.350.000,00, o que 


remete a um retorno de investimento da ordem de 2,1 anos. 


41 CONCLUSÃO E DISCUSSÃO 


Os resultados apresentados pela metodologia mostram que esta tese está correta, 
ou seja, o aproveitamento do calor residual oriundo do condensador de um chiller da central 
de água gelada de um hospital para a geração de água quente é vantajoso em termos 
econômicos e de sustentabilidade. Os resultados da metodologia mostraram uma redução 
de custos operacionais com insumos na central de utilidades de 42%. Por fim, foi feito 
o cálculo da redução de emissão de CO2 com o uso de um chiller tipo heat pump. Uma 
avaliação da mudança de matriz energética de produção de água quente indica que a 
recuperação residual produz uma redução da emissão anual de toneladas de CO2 de 79%. 

Também ocorre uma forte redução do consumo de água da torre de resfriamento, 
da ordem de 55%, que gera uma economia com tratamento químico da água da torre, 
simultaneamente diminuindo a demanda térmica de frio e calor e reduzindo o uso dos 
queimadores de combustível fóssil, além de trazer todos os demais benefícios citados 
acima, tornando esse processo mais sustentável. Ainda conclui-se que quanto maior esta 
simultaneidade, maior a redução do uso de combustível fóssil. 

Foi possível comprovar a viabilidade técnica do uso desse sistema quando existe 
uma simultaneidade de consumo de água gelada e água quente, muito comum em edifícios 
de uso hospitalar. Por fim, deve-se esclarecer que a solução em nada prejudica a geração 
de água gelada, pois o resfriador de líquido terá como função principal a geração de água 
quente, permanecendo o chiller atual com a função de geração de água gelada e setpoint 
independente. 

O relatório Cenário dos Hospitais no Brasil em 2018 (FBH: CNSAÚDE, 2018), 
da Federação Brasileira dos Hospitais, indica que existem cerca de 90.000 leitos que 
se utilizam de centrais de água gelada como sistemas de resfriamento do ar, com uma 
ocupação média de leitos de 70% e um consumo médio de água quente de 125l/dia. 

Assim, extrapolando os Índices resultantes de nosso estudo para o cenário 
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apresentado acima, o presente artigo servirá de base para outros estudos de redução de 
custos operacionais em hospitais e indústrias que são consumidores intensivos de água 
gelada e água quente em processos térmicos, e que atualmente simplesmente descarregam 
potencial térmico para a atmosfera e deixam de reduzir a queima de combustível fóssil e o 
consumo de água. 
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RESUMO: As microalgas são muito utilizadas em 
várias áreas, como, por exemplo, nas indústrias 
alimentar, química e farmacêutica, dentre outras. 
A recuperação de microalgas por meio de um 
processo de floculação-sedi entação seguida 
da microfiltraç o e centrifugação é julgada 
dispendiosa, visto que, principalmente na 
microfiltração, são utilizados materiais onerosos, 
e as membranas precisam ser repostas devido 
ao entupimento que ocorre durante o processo. 
Com intuito de minimizar os valores para realizar 
esse processo, empregou-se uma membrana 
de baixo custo, e assim, usando equações da 
literatura, foi estimado um custo total para o 
módulo de microfiltração, incluindo gastos com 
tubos e válvulas, instrumentos e controles, 
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tanques e estruturas, bomba e membrana. Visto 
que os números encontrados correspondem à 
escala de bancada, foi feito um escalonamento 
para conhecê-los na escala piloto. O custo total 
foi de 1149 dólares para a escala de bancada e 
de 7173 dólares para a escala piloto, o que torna 
o processo viável, diante do fato que no módulo 
microfiltração a estrutura tem uma longa vida útil 
e apenas a membrana seria trocada ao longo dos 
usos. 

PALAVRAS-CHAVE:  Microfiltração, 
escala de bancada e piloto. 


custos, 


ESTIMATED CAPITAL COST OF 
MICROFILTRATION IN À CONTINUOUS 
MICROALGAE RECOVERY PROCESS 

ABSTRACT: Microalgae are widely used 
in various areas, for example in the food, 
chemical, and pharmaceutical industries, among 
others. The recovery of microalgae through a 
flocculation-sedimentation process followed by 
microfiltration and centrifugation is considered 
expensive since, especially in microfiltration, 
costly materials are used and membranes need 
to be replaced. In order to minimize the values 
to perform this process, a low-cost membrane 
was used, and thus, using equations from the 
literature, a total cost for the microfiltration 
module was estimated, including expenses for 
tubes and valves, instruments and controls, tanks 
and structures, pump and membrane. The first 
cost value presented corresponds to the bench 
scale, and a scaling up was also made to obtain 
the pilot scale values. The total cost was $ 1149 
for the bench scale and $ 7173 for the pilot scale, 
which makes the process feasible, since in the 
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microfiltration module the structure has a long life cycle and only the membrane would be 
changed over the cycles of use. 
KEYWORDS: Microfiltration, costs, bench scale and pilot 


11 INTRODUÇÃO 


As microalgas são microrganismos unicelulares encontrados em vários ambientes, 
não só aquáticos, mas também terrestres e representam uma grande quantidade de 
espécies. É estimado que existem mais de 50.000 espécies, mas apenas um número 
limitado, de cerca de 30.000, foram estudados e analisados (MATA; MARTINS; CAETANO, 
2010). Esse organismo, além de ser muito importante para a manutenção da vida nos 
ecossistemas, já que é responsável pela produção de boa parte de O,, também agrega 
valor na área científica e industrial, por isso tem sido bastante procurado para desempenhar 
aplicações biotecnológicas. 

O rápido crescimento das microalgas unicelulares é uma das principais vantagens 
para que as microalgas sejam fundamentais em vários processos. Pigmentos, ácidos graxos 
insaturados, proteínas, polissacarídeos, antioxidantes, enzimas e biodiesel são alguns 
produtos de alto valor formados por meio de técnicas com as microalgas. Ademais, elas 
também são usadas para a produção de ração para gado, farinha de peixe e fertilizantes 
orgânicos. Portanto percebe-se o real valor desse ser vivo para a sociedade (DELATTRE 
et al., 2016). 

Como o uso das microalgas requer que sua biomassa seja separada da suspensão 
primeiramente, algumas técnicas são usadas para esse fim. Contudo, a recuperação é um 
grande desafio devido ao pequeno tamanho das células e seu desenvolvimento acelerado, 
(HENRIQUES, 2017). O método adequado de colheita para remover grandes quantidades 
de água e resultar em biomassa, pode envolver uma ou mais etapas de separação sólido- 
líquido e contribuir com 20-30% do custo total da separação (AHMAD et al., 2012). 

A sequência floculação-sedimentação, seguida da microfiltração e por último a 
centrifugação é uma forma de drenagem de cultura que viabiliza a utilização das microalgas 
em processos futuros. A etapa de floculação-sedimentação é considerada de baixo custo, 
visto que utiliza materiais baratos e como é preliminar, reduz o consumo de energia e custos 
das fases seguintes. Por sua vez, a centrifugação pode ser aplicada a diferentes tipos de 
microalgas e possui uma recuperação rápida e eficiente das células, porém demanda altos 
custos de investimento e de operação (HENRIQUES, 2017). 

Amicrofiltração é um tipo especial de filtração que utiliza como tamanho dos poros da 
membrana o diâmetro de 0,1-10 um, e não necessita de altas pressões como força motriz, 
uma vez que as membranas são relativamente abertas (HENRIQUES, 2017). Além disso, 
esse tipo de filtração utilizando membranas é de fácil operação, baixo investimento, abrasão 
insignificante, baixo consumo de energia, altas taxas de filtragem e sem necessidade de 
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adição química (MENDES, 2019). O sistema utilizado nessa etapa é chamado cross-flow 
ou fluxo cruzado, que possibilita um processo contínuo e mais eficiente, porém possui 
como dificuldades o mecanismo de fouling e a polarização de concentração, no entanto, 
as vantagens dessa maneira de filtrar vencem as desvantagens e tornam o processo mais 
viável. Por outro lado, a filtração convencional pode causar danos ao organismos e sua 
capacidade de processamento de volume é limitada (AHMAD et al., 2012). 

As membranas de microfiltração são recomendadas na recuperação de microalgas 
pelo motivo de proporcionarem maior permeabilidade e menor quantidade de incrustações 
em relação às outras membranas que requerem altas pressões ou operem pelo modo 
convencional (DAS; CHAKRABARTY; BARKAKATI, 2017; HENRIQUES, 2017). Amembrana 
utilizada na microfiltração é a parte mais importante de todo o processo, porque o sucesso 
da separação depende da qualidade dela, ou seja, dos seus componentes e poros. 

Há vários tipos de membranas que podem ser usadas para realizar essa técnica, 
as quais são feitas com agentes porogênicos, capazes de formar pequenos orifícios 
que facilitam a filtração. Dentre elas está a polimérica que por sua vez é formada por 
materiais orgânicos como fluoreto de polivinilideno (PVDF), acetato de celulose (CA), 
politetrafluoroetileno (PTFE), polipropileno (PP), poli (éter sulfonas) (PES), cloreto de 
polivinil (PVC), poliacrilonitrila (PAN) e tereftalato de polietileno (PET). Outro tipo de 
membrana que pode ser utilizada é aquela formada por materiais inorgânicos (cerâmicos), 
dentre eles destaca-se o Al,O, (óxido de alumínio), a alumina, a sílica, o zircônio e o titânio, 
que formam a matéria prima sintética (BILAD et al., 2014). Essas duas categorias se 
enquadram no grupo das membranas comerciais, as quais são constituídas por materiais 
amplamente encontrados no mercado. 

Os materiais cerâmicos apresentam como características a possibilidade de operar 
em uma grande faixa de temperatura, altas permeabilidade e inércia química e baixas 
condutividade térmica e densidade. Essas vantagens propiciam a utilização desses 
materiais na produção de membranas e em várias outras aplicações (HENRIQUES, 2017). 

Ao comparar as membranas poliméricas com as cerâmicas, nota-se diferença no 
que diz respeito ao valor, que para a segunda é julgado mais alto. Por outro lado, as 
membranas cerâmicas (inorgânicas) permitem limpezas mais eficientes e possuem vida 
útil maior do que as membranas orgânicas, além de poder ser utilizadas em sistemas com 
elevada temperatura e pH extremos (TRINDADE, 2010). 

Há pouca literatura e experiência no sentido de estudo dos custos de se utilizar o 
módulo de microfiltração com o objetivo de recuperar microalgas. O principal motivo disso 
é que esse processo envolve um grande número de variáveis e fatores no projeto e na 
sua operação (SETHI; WIESNER, 2000). Por isso alguns custos foram ignorados neste 
estudo e somente aqueles relevantes foram calculados separadamente. Logo, é de grande 
importância analisar o custo de capital gerado no decorrer do método, porque só assim é 
possível julgar a viabilidade dele. 
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Os custos para recuperar microalgas são considerados altos, já que a microfiltraçã 
geralmente utiliza matérias primas sintéticas como forma de construir uma membrana 
cerâmica, isso incorpora a parte dispendiosa de toda técnica (HENRIQUES, 2017). Com a 
finalidade de minimizar os gastos do processo, os componentes da membrana podem ser 
trocados por materiais de baixo custo, como amido de mandioca e casca de ovo. Dessa 
forma, o custo de capital diminui, tornando o trabalho de decréscimo da concentração de 
microalgas mais acessível. 

Existem três tipos de escala possíveis para avaliar o custo e verificar a viabilidade. 
Dentre eles está a escala de bancada, que também é conhecida como de laboratório, a qual 
utiliza equipamentos mais simples e quantidades menores. Nesse estágio é necessário que 
a previsão de desenvolvimento do processo seja precisa para que os custos possam ser 
avaliados de forma realista futuramente (TENNISON, 2000). Já a escala piloto, que mostra 
com mais clareza a praticabilidade, utiliza quantidades satisfatórias e equipamentos que 
reproduzam aqueles utilizados na escala industrial, que por sua vez utiliza valores grandes 
e reais. Neste trabalho, os custos serão apresentados em escala de bancada e piloto. 

Por fim, pode-se perceber a importância em estudar o custo de capital do módulo 
de microfiltração, já que os produtos do processo (microalgas recuperadas) são utilizados 
em várias áreas da química e gestão ambiental, como, por exemplo, no tratamento de 
efluentes 

A microfiltração faz parte da sequência de recuperação de microalgas: floculação 
sedimentação, microfiltração e centrifugação. Sendo assim, a análise do custo de capital 


no módulo de microfiltração é o objetivo principal deste trabalho 


21 MATERIAL E MÉTODOS 


2.1 Material 


As microalgas utilizadas na microfiltração tiveram seu cultivo desenvolvido pelo 
grupo de pesquisa NPDEAS (Núcleo de Pesquisa e Desenvolvimento de Energias 
Autossustentáveis) na Universidade Federal do Paraná. Foi utilizado pelo grupo o meio de 
cultura sintético Chu modificado e. Ademais, utilizou-se como concentração das microalgas 
o valor de 60 mg/L (MENDES, 2019). Os custos da microfiltração no processo de floculação 
sedimentação/microfiltração centrifugação foram avaliados a partir de dados do trabalho de 
MENDES, (2019). 

As características da membrana e os equipamentos utilizados na unidade 
experimental foram retiradas dessa pesquisa. parte de equações foi feita a partir do 
trabalho de (SETHI; WIESNER, 2000) e de (NANDI; UPPALURI; PURKAIT, 2009). 

A filtração utilizada no processo é do tipo cross-flow ou fluxo cruzado e consiste em 
alimentar tangencialmente a superfície da membrana permitindo a permeação, ou seja, 
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o fluxo alimentado passa no processo e se divide em retido e permeado como mostrado 
na Figura 1 item a. Como os poros da membrana vão sendo bloqueados com o passar do 


tempo o fluxo diminui como pode ser observado na Figura 1 item 


Alimentação Retido 4 
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Figura 1 — Filtração tangencial. 


Fonte: (LABVIRTUAL, 2020). 


(HENRIQUES, 2017). 


Tempo 


O sistema de filtração, especificamente, consiste de uma bomba para circulação do 
fluido (marca SEAFLO, modelo SFDP1-014-100-22, vazão volumétrica e pressão máxima 
de 5,0 L/min e 6,9 bar, respectivamente), acionada por fonte (ICEL Manaus modelo PS- 
6000 30V/6A), rotâmetro (Tecnofluid série TRP-330-5), manômetro (963 Glicerina Genebre 
- 3822N com escalas de O a 2,0 bar) e um módulo metálico, chamado módulo de filtração 


onde a membrana cerâmica é posicionada para a realização dos testes (MENDES, 2019). 
A Figura 2 ilustra o sistema de microfiltração que foi utilizado na recuperação das 


microalgas baseado no trabalho de MENDES, (2019). 
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Figura 2 — Sistema de microfiltração 


Fonte: MENDES, (2019). 


A membrana utilizada a partir do estudo de MENDES, (2019) é formada por argila 
pré-tratada a 500 ºC (95% p/p), contendo amido catiônico (2,5% p/p) e casca de ovo (2,5% 
p/p), a qual foi denominada de M1 no referido trabalho. Suas dimensões são apresentadas 
na Tabela 1 e sua forma na Figura 3. A área interna da membrana (A, ..,) para a escala de 
bancada foi calculada pela equação 


Amen = Tt del (1 


na qual, TT é número pi, d é o diâmetro interno e | é o comprimento da membrana. 


Comprimento (m) Diâmetro interno (m) | Diâmetro externo (m) Área interna (m?) 
0,175 0,01075 0,01885 0,005910 


Tabela 1: Dimensões da membrana utilizada. 


Fonte: MENDES (2019). 
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Figura 3: Membrana utilizada por MENDES, (2019). 
Fonte: MENDES, (2019). 


2.2 Métodos 


Com a finalidade de calcular o custo total do módulo de microfiltração tanto para 
a escala de bancada quanto para a escala piloto foram utilizadas equações que geram 
os custos de tubos e válvulas, instrumentos e controles, tanques e estruturas, bomba e 
membrana, separadamente. Após realizar os cálculos com essas equações, foi feita uma 
soma para obter o valor total do processo. A Tabela 2 mostra as equações usadas para 
determinar os custos parciais e o custo total. 


Custos Equações Referências 
(SETHI; 
E E C=h+ (Ane )O (2) WIESNER, 
2000) 
à - (SETHI; 
Custo de tubos e Cry = 5926,13 + (A Po (3) WIESNER 
válvulas ' 2000) - 
Custo de a . 0.66 (SETHI; 
instrumentos e Li = 1445,50 * (Amem) a (4) WIESNER, 
controles 2000) 
SETHI; 
Custo de tanques O ET ; 5% ( : 
 ssias Cre = 3047,21 + (Amem) (5) nas 
E (SETHI; 
Custo da bomba Câombo = tefefae Lo BLETo (Do Py (6) WIESNER, 
2000) 
Custo anual AE a efe (i++ jiu (NANDI; 
do médulo da Coen = ttmeml O * Quem "LO (AIN 1 | uppaLUR 
membrana (+ lu =] PURKAIT, 2009) 
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Custo total O ; : E » ' 
do módulo de Listas E Ery + toc dt Era t Ep t Lhomba + Emem | (8) 
microfiltraçã 


Tabela 2: Equações usadas para determinar os custos parciais e o custo total. 


A partir das equações definidas na Tabela 2, a qual inclui o custo total que diz 
respeito a soma de todos os outros custos, alguns parâmetros, explícitos na Tabela 3, cujo 
valor não muda, foram estabelecidos e utilizados nessas equações para obter resultados 
precisos dos custos. Além disso, é necessário saber outros valores que variam de acordo 
com a escala desejada para calcular os custos. 


Na escala de bancada foram utilizados para A (área interna da membrana em 


m?), Q (vazão em L/min), P (diferencial de pressão em bar), os seguintes valores: 0,005910 
m2, 0,6 L/min e 0,7 bar, respectivamente. 
Já para a escala piloto foi feito um novo cálculo para descobrir o valor da área 


interna da membrana fazendo um aumento de escala e utilizando a seguinte equação: 


td td 
Amemp ES 
E 


na qual Q é a vazão em mº/dia e J é o fluxo no estado estacionário em kg/m2h. Para 
essas variáveis foram usados os seguintes valores: 1 mº/dia, 120 kg/m?h, respectivamente. 
O resultado para a área foi aproximadamente 0,35 m?. A fim de realizar esse cálculo foi 
necessário realizar conversões de unidades e, considerou-se a densidade da água (1000 
kg/m?) como sendo aproximadamente igual à densidade da solução de microalgas, já que 
essa é muito diluída e que o processo atua em regime de estado estacionário (MENDES, 
2019). 

Com o conhecimento desses parâmetros e dessas variáveis, foi possível resolver 


as equações. 


Parâmetro Símbolo Valor 
Índice de atualização do custo 1 2,28 
Fator de ajuste para construção da 
f 1,4 
bomba 1 
Fator a ser ajustado para faixa de f 1 
pressão 2 
Fator correspondente à custos de L 1 
laboratório 
Fator correspondente ao custo do Q 400 
módulo da membrana ($/m?) mem 
Taxa de juros anual Í 0,1 
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Vida útil (anos) a 0,25 


M 


Tabela 3: Parâmetros das equações. 


Fonte: (NANDI; UPPALURI; PURKAIT, 2009; SETHI; WIESNER, 2000). 


O modelo de equações e parâmetros apresentado acima refere-se aos custos 
da operação unitária microfiltração. Com o intuito de estimá-los, o programa Excel foi 
usado para realizar os cálculos, já que por meio dessa ferramenta, os valores obtidos se 
aproximam melhor da exatidão. 


31 RESULTADOS E DISCUSSÃO 


3.1 Escala de bancada 


A primeira condição na qual os custos foram avaliados é a escala de bancada, que 
consiste principalmente em obter o valor do módulo de microfiltração utilizando apenas 
uma membrana com dimensões específicas e usando dados levantados no laboratório. 
Essa forma de realizar os cálculos ainda que em pequenas quantidades, possibilita adquirir 
noções do processo para que sua viabilidade seja julgada. 

Após a estipulação das equações e de seus parâmetros e variáveis, é possível 
realizar o cálculo dos custos. O Excel foi a ferramenta utilizada para esse fim e por meio 
dele obteve-se os valores apresentados na Tabela 4. 


Tipo de custo Valor do custo em dólar ($) 
Custo de tubos e válvulas 687 
Custo de instrumentos e controles 49 
Custo de tanques e estruturas 201 
Custo da bomba 185 
Custo anual do módulo da membrana 28 
Custo total do módulo de microfiltraçã 1149 


Tabela 4: Resultados em escala de bancada (custos parciais e total) 


Fonte: Autora, 2020. 


A partir dos resultados apresentados na Tabela 4, pode-se observar que o menor 
custo foi o da membrana. Isso deve-se ao fato de que os seus componentes são baratos, 
como a casca de ovo e o amido de mandioca. Por outro lado, tudo que envolve a operação 
da microfiltração mostrou-se como custos mais altos, como os tubos, válvulas, tanques 
e estruturas. O motivo disso é a utilização de materiais e equipamentos mais caros para 
construir o módulo de filtração. Embora deva se considerar que a vida útil do módulo de 


filtração é muito longa, em torno de vários anos. E que nesse sistema de microfiltraçã 
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apenas a membrana seria reposta com maior frequência a qual apresenta um custo mais 
baixo. Logo, os resultados possuem uma diferença de valores de acordo com cada tipo de 
custo. 


3.2 Escala piloto 


A partir da dedução de que o processo em questão se comportaria da mesma forma 
para uma escala maior, foram definidos os resultados para a escala piloto, na qual os dados 
levantados no laboratório foram submetidos a um aumento de escala (escalonamento). 
Dessa forma é ainda mais fácil julgar a viabilidade do processo. 

Assim como para a escala de bancada, a ferramenta Excel foi utilizada para o 
cálculo dos custos na escala piloto, que são apresentados na Tabela 5. 


Tipo de custo Valor do custo em dólar ($) 
Custo de tubos e válvulas 3800 
Custo de instrumentos e controles 719 
Custo de tanques e estruturas 1739 
Custo da bomba 185 
Custo anual do módulo da membrana 729 
Custo total do módulo de microfiltraçã 7173 


Tabela 5: Resultados em piloto (custos parciais e total) 


Fonte: Autora, 2020. 


Da mesma maneira que foi observado para a escala de bancada, o custo da 
membrana é considerado bem mais baixo do que àqueles correspondentes à operação 
do processo, devido aos materiais baratos utilizados para sua confecção. Isso ocorre 
porque o que muda de uma escala para a outra é a quantidade de membranas e não sua 
composição. Ao dividir 0,35 m? (área das membranas em escala piloto) por 0,005910 m? 
(área de uma única membrana, utilizada em escala de bancada), o valor retornado é de 
aproximadamente 59, ou seja, foram utilizadas 59 membranas em escala piloto para uma 
vazão de 1 mí/dia. 

Ao comparar os dois tipos de escala no que diz respeito à vazão Q utilizada, percebe- 
se que o escalonamento realizado pode ser eficiente, visto que em escala de bancada a 
vazão foi de 0,6 L/min, que ao converter as unidades, resulta em 0,864 mí/dia. Esse valor 
está próximo de 1 mídia, o qual foi adotado na escala piloto. 

Além disso, outro ponto possível de análise é o custo em dólar gasta por metro 
quadrado de área da membrana, que no caso é de aproximadamente 4738 $/m?, já que o 
valor encontrado de 28 dólares para a escala de bancada é correspondente à 0,005910 m?. 
Com esse dado em mãos é possível saber a quantidade de membranas gastas para 1 m? 


de área da membrana. Ao dividir esse valor por 0,005910 m2 o resultado é 169, ou seja, 1 


Engenharia química: Desenvolvimento de novos processos e produtos 2 Capítulo 7 LES 


m? rende 169 membranas, considerando 0,005910 m? a área de cada membrana. 

Como são muitas membranas, uma grande quantidade de solução pode ser 
filtrada. Como 59 membranas filtram 1 mº/dia de solução, então 169 membranas filtra 
aproximadamente 3 mí/dia. Percebe-se que o valor para a vazão tratada diariamente é 
bastante alto, por isso seria viável utilizar a composição da membrana adotada (argila 
pré-tratada, amido catiônico e casca de ovo), já que ela é de baixo custo por ser feita com 
resíduos. 

Na literatura existem poucos estudos relacionados ao custo de operação do 
processo em questão. No entanto, no que diz respeito ao valor para a membrana, incluindo 
a fabricação e os materiais utilizados, alguns autores apresentaram seus resultados. 

Estudos realizados por Nandi et al. (2009) mostraram 130 $/m? como valor da 
membrana de baixo custo. Esse valor encontra-se abaixo da média do custo de uma 
membrana desenvolvida com materiais cerâmicos, que segundo seu trabalho, é de 400 $/ 
m2. Outros autores como Tennison (2000) e Das et al. (2017) chegaram em seus estudos, 
para os custos da membrana, nos valores de 598,90 $/m? e 92 $/m?, respectivamente, 
sendo o último desenvolvido para uma membrana feita com materiais crus, ou seja, 
materiais de baixo custo, assim como foi feito nesse trabalho. 

Como pode ser visto, o valor encontrado por metro quadrado (4738 $/m?) é bem 
acima dos valores dos trabalhos citados, embora espere-se que com o aumento de escala, 
os custos de produção das membranas diminuam. No entanto, vale ressaltar que nesses 
estudos, não é apresentado o comportamento das membranas em relação ao fluxo, isto 
é, a eficiência das membranas. Por isso torna-se inviável comparar esses custos de forma 
realista. 


41 CONCLUSÃO 


A partir dos resultados obtidos, pode-se concluir que a microfiltração utilizando 
membranas de baixo custo na sequência de recuperação de microalgas: floculação 
sedimentação, microfiltração e centrifugação é viável. Isso pode ser observado por meio da 
análise do custo de capital para esse processo, que gerou um custo total de 1149 dólares 
para a escala de bancada e de 7173 dólares para a escala piloto, ao somar todos os custos 
parciais. 

As equações adotadas para a análise do custo de capital foram selecionadas 
conforme o tipo de custo. Elas foram divididas em custos de tubos e válvulas, instrumentos 
e controles, tanques e estruturas, bomba, membrana e a equação de custo total, que é 
a soma dos outros valores. A partir disso, foi possível determinar o resultado para cada 
equação, utilizando a área interna da membrana, calculada de acordo com a escala. 

Ao comparar os dois tipos de escala, observou-se uma eficiência no escalonamento 


adotado já que os valores para vazão ficaram próximos, os quais foram 0,864 mí/dia e 1 
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mídia, para a escala de bancada e piloto, respectivamente. Além disso, foi observado a 
quantidade de membranas utilizadas e a porção de solução que elas filtram, mostrando que 
o processo de microfiltração para a separação de microalgas é apropriado, porque o valor 
que é filtrado é alto 

Referências na literatura são escassas para um custo total, mas não para a 
membrana. Como os estudos levados em questão não apresentam o comportamento das 
membranas, os valores encontrados não foram apropriados para uma comparação mais 
realista. 

Portanto, a análise do custo de capital possibilitou que a etapa de microfiltração no 


processo tratado fosse factível e possível de ser aplicada. 
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RESUMEN: Se llevó a cabo la remoción de 
arsénico (3+), por medio de la adsorción sobre un 
sustrato proveniente de la pirolisis de Ilantas fura 
de uso (carbón de pirolisis), debido al contenido 
de azufre y otros metales, fue necesario 
realizar su acondicionamiento por medio de una 
desulfuración oxidativa (ODS) y un tratamiento 
químico. Las concentraciones de arsénico 
fueron medidas mediante espectrofotometria 
UV/Vis antes y después de cada tratamiento 
por medio de espectrofotometria UV/Vis; para 
comparar el rendimiento del carbón de pirolisis 
se llevaron pruebas similares con un carbón 
comercial de las mismas características. Se 
realizaron curvas de equilibrio y se propusieron 
dos modelos de adsorción para realizar el ajuste 
de los datos experimentales. Se logró disminuir 
considerablemente el azufre presente en el 
carbón mediante la ODS, sin embargo no se 
logró su eliminación, lo cual sería una limitante 
para la aplicación del carbón en el tratamiento del 
agua. El carbón de pirolisis tuvo un rendimiento 
superior al carbón comercial, se encontró un 
modelo tipo Lagmiur para los datos de equilibrio, 
la adsorción de arsénico sobre el carbón de 
pirolisis resulto ser favorable y representa un 
medio sustentable para la remoción de metales 
altamente contaminantes. 


INTRODUCCIÓN 

Debido a su alta toxicidad, no bio- 
degradabilidad y acumulación a través de la 
cadena alimenticia, el arsénico y otros metales 
pesados resultan ser muy problemáticos y 
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amenazantes para el medio ambiente y los seres humanos. Debido a esto, la Organización 
Mundial de la Salud, recomendó que el máximo permisible debe ser de 10 ug/L para el 
arsénico y para cromo el 0.05 mg/L en agua para consumo humano. Por lo tanto es muy 
importante lidiar con las aguas residuales contaminadas antes de su descarga al suelo y 
a los cuerpos de aguas. Muchos investigadores han empleado la precipitación química, 
separación por gravedad, ósmosis inversa, intercambio iónico, electrocoagulación, 
electrodiálisis, flotación y adsorción para la remoción de fosforo y metales pesados de la 
red del alcantarillado [1]. Recientemente el negro de carbono proveniente de la gasificació 

del aceite crudo de los fondos de refinería, fue acondicionado para adsorber colorantes, 
mostrando una gran capacidad de adsorción [2]. Los materiales de ingeniería resultan 
ser una fuente alterna para el tratamiento de metales traza, de ahí la importancia de este 


trabajo para elaborar un material a partir de la pirolisis de llantas post-consumo [3,4,5]. 


METODOLOGÍA 


Reactivos y Materiales. Se empleó arsénico (As,O,) de la marca Sigma aldrich 
(CAS 1325-53-3) el cual fue empleado tal como fue recibido. 

Preparación del carbón del pirolisis. El carbón de pirolisis fue desulfurado mediante 
peróxido de hidrogeno (10%) y desmineralizado mediante ácido clorhídrico 1 M. El carbón 
desulfurado fue lavado con metanol para remover los compuestos oxidados y posteriormente 
lavado con suficiente agua desionizada. Inicialmente el carbón fue colocado en una solución 
de HCI por un tiempo de 30 minutos en agitación constante, posteriormente fue lavado con 
abundante agua y secado a 80 “C, para posteriormente Ilevar a cabo la desulfuración. 

Caracterización del carbón. El carbón fue caracterizado mediante patrones de 
difracción de rayos X, análisis químico elemental (fluoresencia de rayos X) y área superficia 
(BET), antes y después del tratamiento de desulfuración y del tratamiento químico con HCl. 

Análisis de arsénico (3+)6]. Para el análisis del As (3+) se empleó el método 
espectrofotométrico con dietil ditio carbamato de plata, en este método el arsénico presente 
se reduce a arsina por medio de zinc metálico en solución acida, la arsina formada reacciona 
con el dietil ditio carbamato de plata para formar un complejo rojo soluble, intercambio 
del As por la plata (535 nm), cuya coloración es proporcional al arsénico presente en la 
muestra (reacciones 1 y 2). Se utilizó un espectrofotómetro Varian Cary 50 C UV/Vis para 
determinar las concentraciones del complejo rojo al inicio y al final de cada tratamiento 
ZA Zn + OH! = Jncl, + 2Azh, Ro, 

AsH, + AgSCSN (CM), — Às — Complejo (coloración roja) | 

Desemperfio del carbón de pirolisis. Para la comparación de resultados, se utilizó 
carbón activado comercial GAC 1240W, para procesamiento de agua potable, activado por 
vapor granular marca Norit'M (GAC). 


Curvas de equilibrio. Se realizaron curvas de equilibrio, el tiempo establecido fue 
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de 180 min para alcanzar las condiciones de equilibrio, las concentraciones en el equilibrio 


fueron determinadas mediante la ecuación 1 y se estableció un modelo de acuerdo a la 


ecuación 2. 
Ci— 
tp = cd La (o) 
me 
RR (é) 
Ing, = Ink 4 BA, Ines (3) 
RESULTADOS 


Análisis químico Elemental. En la figura 1 y 2 se muestran los resultados del análisis 
químico elemental antes y después del tratamiento. Los elementos que se presentaron 
en mayor proporción fueron zinc, azufre y hierro, no se encontró la presencia de vanadio, 
cadmio ni otro elemento traza, aunque si se detectó aluminio al inicio del tratamiento. 
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Figura 1.0 Análisis químico elemental después Figura 2.0 Comparación de los tratamientos en la 
del tratamiento químico y desulfuración. remoción de los elementos. 


Área Superficial. En la tabla 1 se muestran los resultados de la fisisorción con 
nitrógeno gaseoso sobre el carbón de pirolisis antes y después de los tratamientos, debido 
al diámetro del poro el material se puede clasificar como un material meso poroso, lo cual 
refleja espacios disponibles para la adsorción. La desulfuración y la desmineralización 


aumentan el área superficial y disminuyen el volumen del poro 


Carbón Carbón desmineralizado 
desmineralizado y desulfurado 
Área superficial ( 2/9) 51.77 63.07 
Volumen del poro (cc/g) 0.202 0.173 
Diámetro del poro (À) 14.288 14.288 


Tabla 1.0 Caracterización del carbón de pirolisis mediante análisis BET. 
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Difracción de rayos X. En la figura 3 se muestra el patrón de difracción en el carbón 
antes y después del tratamiento, solo se observó respuesta para el Zn (cristalino) que 
puede ser en su forma más común (sfarelita hexagonal). 
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Figura 4. Efecto de la concentración inicial en 
la remoción de As (+3) en carbón de pirolisis 
(W=0.05 9) 
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Figura 3.0 Patrones de difracción de rayos X Figura 5. Efecto del tipo de carbón (1: carbón de 
antes y después del tratamiento. pirolisis, 2: Carbón comercial Norit TM) (Co=450 


mg/L, W=0.05 9) 


Pruebas de adsorción con el carbón comercial y el de pirolisis. En la figura 4 se puede 
observar el efecto de la concentración inicial de As(3+) sobre la remoción en el carbón 
de pirolisis, la concentración inicial afecta fuertemente la remoción, siendo mayor cuando 
se reduce la concentración inicial (1000<820<635<450 mg/L), a partir de los 15 minutos 
de contacto con el carbón se alcanzó un 81% de remoción para la menor concentración 
inicial estudiada, a partir de este tiempo la remoción permanece constante. En la figura 5 
se muestra el efecto del tipo de carbón en la remoción de As(3+), el carbón 1 (carbón de 
pirolisis) tuvo un desempefio superior al comercial (2), para un tiempo de 15 minutos este 
alcanzó únicamente un 33.5% de remoción en comparación con el carbón de pirolisis que 
alcanzó un 81.0% para este mismo tiempo. 
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CONCLUSIONES 


El carbón de pirolisis representa una opción para el tratamiento sustentable del 
agua contaminada con elementos tóxicos como el arsénico, sin embargo este puede 
verse limitado por la presencia de otros elementos tóxicos como el azufre contenido en 
el carbón de pirolisis por lo que se tienen que realizar la ultra-desulfuración (no realizada 
en este estudio), para asegurar su probable aplicación. El modelo de primer orden para 
las concentraciones en el equilibrio (n=0.996) representó adecuadamente los datos 


experimentales y aseguró una adsorción favorable (k=45 mg/mg). 
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